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AN ANALYTICAL SATELLITE ORBIT PREDICTOR (ASOP) 
By Stephen E. Starke 

Analytical and Computational Mathematics, Inc. 


1,0 INTRODUCTION 

This report contains the documentation and user’s guide for the Analytical 
Satellite Orbit Predictor (ASOP) computer program. ASOP is based on mathema- 
tical methods that represent a new state-of-the-art for rapid orbit computa- 
tion techniques. The theoretical development of these methods has been 
carried out during the past few years, and they are now being put into the 
form of an operational computer program. ASOP is intended to be used for 
computation of near-Earth orbits including those of the Shuttle/Orbiter and 
its payloads. 

Orbit computation met lods can usually be given one of the following two 
classifications: 

a. Numerical methods - The calculations are carried out in a step-by-step 
manner. High precision is possible, but computer runtime can be excessive. 

b. Analytical methods - The calculations are carried out in one step regard- 
less of the prediction interval. Therefore, these methods have extrerely 
fast computation times. 

In the past, analytical methods have not been widely used because they were 
less accurate and required much more computer codes than numerical methods. 

The Poincare-Sirai lar Element Method (PS<j)) used in ASOP overcomes these disad- 
vantages. Tt is possible to compute near-Earth orbits bo within an accuracy 
of a few meters. Recursive equations are used instead of complicated formulas. 
Execution time of ASOP is on the order of a few milliseconds, tt is planned 
to iiiclude in ASOP the perturbations due to atmospheric drag, Sun and Moon gra- 
vity and all terms of the geopotential. 

Tn reference 1, Mueller describes the relationship between the PS4) elements and 
the C3rte.sian coordinates and establishes the PS(J) perturbed equations of motion 
These elements and the a.nsociated equations of motion were used by Bond (ref, 

1) to develop a completely analytical solution to the J 2 per’turbablon problem 

Ret enlly, the analytical .solution has been expanded by Schelfele (ref. 3) to 
inr lude the perturbations due to atmospheric drag. Mueller (ref. 4) has devel- 
oped the higher order zonal geopotential terms. Time dependent (tesseral) geo- 
potential terms for the solution are under development. Later, it is expect. ed 
th .t the third body perturbation.s ui‘ the .Sun and Moon will also be added to the 
sol iitlon . 



The current version of A30P Includes the J2 geopocentlal terms. The coding 
is documented in this report. It is planned that the additional perturbations 
mentioned above will be added to ASOP and will be documented in separate reports. 

This ASOP program has been designed in two versions: (1) a stand-alone version 

that can be used interactively to obtain immediate results to a specific prob- 
lem and (2) a user-subroutine package that can be used by other software sys- 
tems. Both versions were designed to be small and to execute quickly. The 
average execution time is ilO milliseconds for the stand-alone program and 5 
milliseconds for the user-subroutine package, while using 3111 10 and 2001io* 
words of storage, respectively. 

Both versions of the ASOP program have been written in UNIVAC standard FORTRAN-V 
and are available to the public on file NUMEG under the qualifier FM6-NO8569 
on the UNIVAC 1110-Exec 8 system. This document is intended to instruct the 
user in the operation of the ASOP program on this machine and to document the 
individual ASOP subroutines. 
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2.0 PS ELEMENT FORMULATION 


The PS (Poincare-Similar) element formulation is deaoribed in the following 
subsections. An exact development of this element set is describee in reference 
2 and a description of the variables can be found in Appendix E. 


2.1 TRANSFORMATION FROM CARTESIAN COORDINATES TO PS ELEMENTS (XTOPS) 

The transformation into the PS element set (a,"?) from a given set of Cartesian 
coordinates (X,V,r) is accomplished within the subroutine XTOPS. It is 
assumed that these coordinates are the initial conditions so that the physical 
time t and the independent variable t both equal zero. The PS elements 
are then computed using the equations 


(T, = tan 


-1 


'Xjj + R"a3 

X7TW3 


Og = Zg cos Oj - sin 


O3 . 0^) 


Oi, = t - 


(2pi, )3/2 


ZgQ 


^E - 41 

Vi - e ^ ^ 


The other required relations are; 


p = 0 - Vg^ + 2r^V + — t- 

pg = Zg sin cos a.j 

P3 = 2Gy/y2(G + G^) 



G=XV-XV , G=XV- XV , G„=XV-XVj 
X y z z y ’ y zx xz’ z xy y^ 


? ? P 

G=.,/G +G +G 
X y z 


p = 1/11 (G - p^ + , q = G - 1/2 Pi + U/2>y2P4 


R* = rR/2G , H = 2GX^/r 


Z, = (p/p - 1)/Q , Zg = rp/(q(2q - G)) 
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Q = 1/M (Pi, 


X'V I 2 I 

r = -— , Vl - = P ^2 Pj,/M 

E - 4» = -2 tan“‘'(Z2Q/(1 + Vn~8 + Z^Q)) 

This transformation is performed only once for a given set of Cartesian 
coordinates . 


2.2 TRANSFORMATION FROM PS ELEMENTS TO CARTESIAN COORDINATES (PSTOX) 

This transformation is performed when any intermediate printout is desired 
or when the final condition is met. Therefore, this transformation should 
be coded with emphasis on the speed of calculation. 

The Cartesian coordinates defined in terms of the PS elements are given by 
the equations; 


X„ = r cos a, - R a_, , V = r cos a. - - sin a 


1 "'3 


'1 


X^ = r sin - R 


V = r sin a, - - cos a. 
y 1 r ^ 


X = R* yj2(G + H) , V = R*V2(G + H) 


The other necessary relationships are 


R*a 


3 


3 


r = p/(l + e cos (j>) I R* = rR/2G 

r = e sin (J)/p(2q - + 1/2 + a^^)] 

G 

R* = (Rr + RrO/2G , R = — (p.^ cos sin 

r^ 

2 

R = sin cos , G = - 1/2 ( p^ 


) 




<1 



H t G - 1/2 , q = - 1/2 

p X 1/y [W\f^~ 1/2 (P2^ + o^^)) 

e sin (|i = Q(p 2 sin cos ) 

e cos (|) B Q(p2 008 - Og sin ) 

Q = 1/y I p^^[2y/Y^- 1/2 ip/ + a/)]] 

The physical time t is computed using 

t B 0. + (E - (|) - - e sin (ji ^/l - e^) 

^ (^P^)3/^ P '' 


p^ - y/yi^) 


(2.2a) 


where an expression for E ■■ s}; is given by 


E - (f = -2 tan**^ |^e sin <{)/(1 + ^1 - e^ + e cos 4i)J 
2.3 TIME TERMINATION PROCEDURE 

Because the PS element set uses the true anomaly t as the independent vari- 
able, an iteration procedure is necessary to stop at a specific time tfinal* 
Within the ASOP program, this iteration is performed by the TIMERS subroutine 
in the following manner: 

An expression for the derivative of time with respect to the true anomaly T 
is given in the PS theory as 


dt 

di 



/q 


This expression can be linearly approximated by the expression 


pt 2 

— = r /q (2.3a) 

At 


ORIGINAL P \GE IS 
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where At « - tj, and At * - Tj,. Equation 2.3a then yields a re- 

cursive formula for refining an initial estimate of t of 

Using equation 2.3b, an initial estimate of t is refined until the associated 
value of t equals the desired final time 

In order to start the iteration, a suitable initial value of x is necessary. 
This value is determined by first assuming that a circular orbit is being 
used. With this assumption, equation 2.2a reduces to 


= t 


and an initial approximation for t can be written as 


3crn 


The value of 30[|/3 t is taxer. ; ron the analytical theory (see reference 2) 
during the initialization '.rnctouro. 

Therefore, the full algorithm, is 

a. Set the iteration counter n to zero, and compute the initial approximation 
for In = To- 

b. Determine the PS elements at r^j. 

c. Determine the time t^ at Tj^* 

d. If - t^l £ TOLerenoe, then STOP; otherwise n = n + 1 

e. If n > n[nax> then print a diagnostic message and STOP 

f. Compute a new approximation for Xf^ using 

Tn = ^n-1 - “ '■final 

g. Go to step b. 
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Values- for TOLerence and n have been preset within the TIMERS subroutine 
-7 max 

to 10 and 5, respectively. 
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3.0 DESCRIPTION AHD STRUCTURE OF A30P 


TL-: A30P program was designed as a top-down structured program consisting of 
11 subroutines (modules) and a main (driver) program. Within this set of sub- 
routines, a rubset of seven subroutines comprises the removable ASOP subroutine 
package. The package has been designed for easy incorporation into existing 
user software. This type of division nas resulted in some redundant operations 
but Joes not seriously affect the overall performance of the ASOP program. 

The two sets of subroutines are described in the following subsections and each 
subroutine is documented in section 3.3. 


3.1 THE ASOP PROGRAM 

Basically, the ASOP program, shown in figure 1, consists cf four segments: 
a main program or driver, an input routine, an output routine, and a removable 
ASOP subroutine package. The removable package is described in section 3-2 
so that on.y the first three segraento will be described here. 

The purpose of the main program is to call all the necessary subroutines (input, 
output, constants, etc.) in a specific sequence to produce the desired results. 
In particular, the main program provides for the repetitive call to the ASOP 
subroutine in order to produce a satellite ephemeris. 

All input to the program is controlled by the subroutine INPUT. Its primary 
functions are 

a. . Set all default values 

b. Accept input from the NAMELIST $INPUT 

c. Convert any input values that require conversion 

d. Issue normal program termination command 

Output from the ASOP Program is performed only by the subroutine OUTPUT. This 
subroutine contains all the FORMAT specifications used for normal output from 
ASOP, The output program will also perform any conversions required to make 
the output more understandable. This involves converting Cartesian coordinates 
to Keplerian elements, radians to degrees, and time units to days. 


^All 


error output is controlled by l.tie 


individual subroutines. 


ORIGINAL PAGE IS 
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3.2 THE ASOP SUBROUTINE PACKAGE 

The ASOP subroutine package has been designed as an independent segment of the 
ASOP program so that it can be incorporated inco existing software with little 
or no modifioations. 

This package, as with the ASOP program, has four basic parts (fig, ?); 

a. A driver subroutine (ASOPj 

b. Coordinate transformations (XTOPS, PSTOX) 

c. A stopping routine (TIMEPS) 

d . The analytical theory (PSANS) 

Along with the five subroutines mentioned, there are two additional subroutines, 
P0TJ2 and CONST, that oerform functions that are used by more than one part so 
that they cannot be specifically Included in one section. 

Subroutine ASOP performs the same task as the main program in that it calls all 
the necessary subroutines as dictated by the input values. In the case of the 
ASOP subroutine, however, the input values are given through an argument list, 
and a few COMMON blocks. Therefore, if a satellite ephemeris is desired, the 
user must supply the necessary coding within his own software. 

The basic input to the subr .Jtine is: 

a. The Cartesian coordinates (X) 

b. A stop value (STOP) 

c. A stop flag (ISTOP) 

d. A new data flag (NEWX) 

This input is fully described in sections 3.3.3 and 11.2.3. 

Because the analytical theory has been developed In PS elements, it is neces- 
sary to perform transformations to and from the element set. Transformation 
into the PS elements from Cartesian coordinates is performed by the subroutine 
XTOPS while the reverse transformation is performed by the subroutine PSTOX. 

A time-stop subroutine has been included because the PS elements use the true 
anomaly as the independent variable and do not use time. This stopping routine 
is an iterative procedure and is described in .section P.3. 

PSANS, the analytical theory subrouiine, is the nucleus of the ASOP program. 

At present, the theory contains only the perturbations due to and these 

equations are given in Appendix F. Hnv:ever, the theory will soon be expanded 
to in' tude other perturbations such as drag (ref. 5), and the higher order geo- 
potential terms. 
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3.3 MODULE DESCRIPTIONS 

This section will give a complete description of the subroutines currently used 
in the ASOP program. Each description will contain a brief statement as to 
the purposv or use, of the subroutine as well as a complete description of all 
variables used within the subroutine. Also included are lists of the named 
COMMON blocks used, external references to other ASOP subroutines, and other 
ASOP subroutines that reference the subroutine being described, information 
is also available as to the calling sequence of the subroutine, subroutine 
size, and the average execution time of the subroutine. Each description is 
followed by a general flow chart of the subroutine (figs, 3 through 14). 

Each program is listed alphabetically, with the exception of the MAIN program, 
which is described first. 


3 . 3.1 MAIN Program (Driver) 

Purpose : Driver for the Analytical Satellite Orbit Predictor (ASOP) program 

Calling sequence : None 

Called by : Operating system 


Subroutlnes/functions called : ASOP, CONST, INPUT, OUTPUT 


Named COMMON : /CARTC / 

/CBASIC/ 
/CPRINT/ 
/END / 
/PS / 


X(8),R 

PI,TWOPI,DEG,RAD,DAy 
PRINT , IPRINT , lELPRT , lUNITS 
ST0P,IST0P 
SIG(8),TAU 


Program data : Size = 124g(84-^Q) words compiled. 

Execution time = 34 milliseconds to 42 milliseconds, depending 
upon the input values for PRINT and IPRINT (does not include 
time for data units conversions and initial condition printout). 


FORTRAN Input/ 

variable Dimension Type output Units Description 


DAY 


IPRINT 


1 DP 


1 r 


sec/day or Conversion of days into hours, 
min/day or minutes, or seconds. Its numer- 
hrs/day leal value is controlled by 

lUNITS. 


I — Flag to determine if the inter- 

mediate printout is to be given 
on days or revolutions. 

= 0 no printout 
= 1 days 
= 2 revolutions 
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FORTRAN Input/ 

variable Dimension Type output Units Desorlption 


ISTOP 

1 

I 

I 


Flag to determine if the pro- 
gram is to terminate after a 
given number of days or 
revolutions as specified by 
STOP. 

= 1 days 
= 2 revolutions 

NEWX 

1 

I 

I 


Flag to determine if the ASOP 
program is to be initialized. 

= 0 no 
= 1 yes 

PRINT 

1 

DP 

I 

days or 
revs 

Increment for which the Inter- 
mediate printout is desired 
(valid only if IPRINT ^ 1) 

STOP 

1 

DP 

I 

days or 
revs 

Value at which the program is 
to stop execution (units are 
determined by ISTOP). 

STOPDT 

1 

DP 

“ 

days or 
revs 

Value at which the next inter- 
mediate printout is desired 
(valid only if IPRINT >1). 

TAU 

1 

DP 

I 

rad 

Independent variable of the PS 
elements; it is defined such 
that REVS = TAU/2TI. 

TWO PI 

1 

DP 

I 

rad 

2tt 

X 

8 

DP 

I 

Defined by 
input val- 
ues (see 
CONST) 

X(1) > X(3) = X 
Xf4) X{6) = t 
X(7) = time 
X(8) = total energy 


ORIGTMAL PAGE IS 
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5T0PDT= PRINT 


CALL ASOP 



YES 


.. , y 

: CALL OUTP, T 
! ; ;rr; ! t P.J’. 


► call ASOP 



FUjure 3.~ MAIN program flow charts 





3.3.? AEIXY2 (Subroutine) 


Purpose : Transform the Keplerian elements (a,e, into Cartesian coordi- 

nates (ti^) 

Calling sequence ; CALL AEIXYZ 
Called by ; INPUT 


Subroutines/funotiona used: None 


Named COMMON ; /CARTC / 
/KEPLER/ 
/KONST / 
/RPOOL / 


X(3),V(3), TIME, ENERGY, R 
A,E,XI,OMEGA,XNODE,XM 
XMU , XMUI , SQTMU , SQTMUI , EPS 

COMEGA , SOMEGA , CNODE , SNODE , SINC , B ( 3 , 3 ) , EA , EAO , SINEA , 
COSEA, X1(3),V1{3) 


Equivalence ; (B{3 , 3) ,CINC) 


Program data : Size = 356g (238 ^q) words compiled 

execution time = < 1 milliseconds 
Subroutine valid only when |e| < 1, 


FORTRAN 

variable 

Dimension 

Type 

Input/ 

output 

Units 

Description 

A 

1 

DP 

I 

see 

appendix A 

Semimajor axis of the orbit 

B 

(3,3) 

DP 

— — 


Keplerian elements to Cartesian 
coordinates transformation 
matrix 

CINC 

1 

DP 

— 


Cosine of the orbital inclina- 
tion with respect to the Earth' 
equator (cos i). 

CNODE 

1 

DP 

— 

— 

Cosine of the argument of the 
ascending node (cos fi). 

COMEGA 

1 

DP 

— 

— 

Cosine of the argument of 
pericenter (cos uJ). 

COSEA 

1 

DP 

— 

— 

Cosine of the eccentric anomaly 
(cos E) . 

E 

1 

DP 

I 

— 

Orbital eccentricity (e). 
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FORTRAN 

variable 

Dimension 

Type 

Input/ 

output 

Units 

Description 

EA 

1 

DP 

— 

rad 

The eccentric anomaly of 
the satellite computed from 
Kepler's equation. 

EAO 

1 

DP 


rad 

Old value of EA; used when 
iterating to solve Kepler's 
equation. 

OMEGA 

1 

DP 

I 

rad 

Argument of pericenter (w). 

R 

1 

DP 

0 

see 

appendix A 

Magnitude of the position 
vector of the satellite. 

SING 

1 

DP 


— 

Sine of the orbital inclina- 
tion with respect to the 
Earth's equator (sin i). 

SINEA 

1 

DP 

— 

— 

Sine of the eccentric anomaly 
(sin E). 

SNODE 

1 

DP 

— 

— 

Sine of the argument of the 
ascending node (sin Ji). 

SOMEGA 

1 

DP 

— 

— 

Sine of the argument of 
pericenter (sinw). 

V 

3 

DP 

0 

see 

appendix A 

Velocity vector of the satel- 
lite with respect to the 
Earth's equator 






V(1) = V 
V(2) = 
V(3) = V^ 

VI 

3 

DP 

— 

see 

appendix A 

Velocity vector with respect 
to the orbital plane 






VK1) = V 

Vl(2) = V'^ 

VK3) = V’^ = 0 
z 

X 

3 

DP 

0 

see 

appendix A 

Position vector of the satellite 
with respect to the Earth's equa' 


X(1) = X 
X{2) = Y 
X(3) = Z 
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FORTRAN 

variable 

Dimension 

Type 

Input/ 

output 

Units 

Description 

XI 

3 

DP 

— 

see 

appendix A 

Position vector of the satellite 
wibh respect to the orbital 
plane 






XK1) = X* 
X1(2) a vt 

XK3) = Z' = 0 

XI 

1 

DP 

I 

rad 

Orbital inclination to the 
Earth's equator (i) 

XM 

1 

DP 

I 

rad 

Mean anomaly of the satellite 
(Ml 

XMU 

1 

DP 

I 

see 

appendix A 

Gravitational constant for the 
central body (p) 

XNODE 

1 

DP 

I 

rad 

Argument of the ascending node 

cn) 


ORIGINAL PAGE IS 
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COMPUTF: 
s 1 n p j c 1. I r- < ' 

CJ J1 6 



COMPUTE: 
TRANCFORMATtON 
matrix B ! 


■ “I’' 

I_. 

SOLVE , 
KEPLER S 
EQUATION 


COMPUTF ; 
r n' V ^ 





- - j 


I 

t 



R f I URN 


/ 


Figure A AEIXYZ flow chart. 
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ASOP (Subroutine) 


Purpose: Driver for the analytical section of the ASOP proprram; It performs 

all the operations required by the analytical program PSAMS 

Calling sequeuce: CALL ASOP (X,STOP,ISTOP,NEWX) 


Called by: MAIN 


Subroutlnes/funotions uaed; CONST, PSANS, PSTOX, TIMEPS, XTOPS 


Named COMMON : /CARTC / 

/CBASIC/ 
/PS / 
/PSTIME/ 


XIN(8),R 

PI.TWOPI,DEG,RAD,DAY 

SIG(8),TAU 

CLO,FAKTPS,XXX 


Program data : Size = 1273 words compiled 

execution time = 2 to 8 tnilHseconda depending on the value of 
ISTOP or the eccentricity of the orbit. 


FORTRAN Input/ 

variable Dimension Type output Units Description 


CLO 

1 

DP 

0 


Initial value of used hv 

the time iteration stopping 
procedure (TIMEPSl. 

DAY 

1 

DP 

I 

see 

appendix A 

Conversion of days into hours, 
minutes, or seconds 

ISTOP 

1 

DP 

I 


Flag to dete>-mine if the value 
of .STOP is given as davs or 
revolutions . 

= 1 days 
= 7 revolutions 

NEWX 

1 

T 

I 


Flag to determine \e ftsop 

program is to be initialized. 

= n no 
= 1 ves 

STG 

8 

DP 



The PS elements 
1 * • "a < ’ ^ll 

Where Pji = ->■ 

STOP 

1 

DP 

I 

days nr 
revs 

Value at which the final state 
vector 1? required; units are 


set hv ISTOP. 


Id 
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fcHTRAN 

variable 

TAU 

TWOPI 

X 


XIN 


Input/ 


Dimension 

Type 

output 

Units 

Description 

1 

DP 

— 

rad 

Independent variable of the PS 
elements. 

1 

DP 

I 

rad 

2ti 

8 

DP 

1/0 

see 

appendix A 

Initiai/final state vector 


X(1) X(3) = X 

X(i4) + X(6) = t 
X(7) = time 
X(8) = total energy 

If initializing NEWX = 1 , 

X(7) and X(8) will be set 
bo 0 

8 DP I/O see Identical to X but allows the 

desription ASOP subroutine to be removed 
of X from the stand alone program. 
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Figure 5.- ASOP flow chart. 
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3-3. CONST (Subroutine) 


Purpose i Initialize the mathematical and physical constants needed to execute 
the ASOP program. 

Calling sequence : CALL CONST 

Called by : ASOP, INPUT 

Subroutines/functions used: None 


Na med COMMON ; /CBASIC/ 
/CPRINT/ 
/GEO / 
/KONST / 
/PERTRB/ 


PI,TWOPI,DEG,RAD,DAY 

X,I(2),IUNITS 

RE,CJ2 

XMC , XMUI , SQTMU , SQTMU I , EPS 
IDRAG.ILONG 


Program data ; Size = 26?8 words compiled 

execution time = < 1 rallliseconda 


FORTRAN 

variable 

Dimension 

Type 

Input/ 

output 

Units 

Description 

CJ2 

1 

DP 

0 

— 

coefficient of the central 

bbdy 

DAY 

1 

DP 

0 

see 

appendix A 

Value to convert days into 
seconds, minutes, or hours 

DAYS 

7 

DP 

— 

see 

appendix A 

Storage array holding the 
possible values of DAY 

DEG 

1 

DP 

0 

deg/rad 

180/71 

EPS 

1 

DP 

0 

— 

3/2 u r/ 

TLONC 

1 

I 

I 


Flag to determine the type 
of geopoLential terms to 
be used. If ILONG = 0 the 
two-body oi'bit is assumed 
(J2 = 0) 

ION ITS 

1 

I 

! 


Flag to determine what units 
are to be unnl for the 
caiculatiins, 

= 1 km, sec 
= 2 nm,sec 
= 3 ft, sec 
= ^ m,sec 
= 'j km, hr 
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FORTH AN 
variable 

Dimension 

Type 

Input/ 

output 

Uni ts 


Description 

lUNITS 

(conc'd) 






= 6 nra,hr 
= 7 E.r. ,min 

PI 

1 

DP 

0 

rad 


71 

RAD 

1 

DP 

0 

rad /deg 


11/180 

RE 

1 

DP 

C 

see 

appendix 

A 

Central body equatorial radius 

RES 

7 

DP 

— 

— 


Array containing the possible 
values of RE 

SQTMU 

1 

DP 

0 

see 

appendix 

A 


SQTMUI 

1 

DP 

0 

see 

appendix 

A 

i/VU' 

TWOPI 

1 

DP 

0 

rad 


2r 

XMU 

1 

DP 

0 

see 

appendix 

A 

Gravitational constant for the 
central body ( p) 

XMUI 

1 

DP 

0 

see 

appendix 

A 

1/H 

XMUS 

7 

DP 

— 

see 

appendix 

A 

Array containing the possible 
values of XMU. 


p-i 
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BEGIN 


<6<|UNlTSi7> 


NO 


1UNITS = 1 


YES 

' 

• W - — . -p 

.. .. .. .. 

---11. -Ill 


SET MATH 

constants; 

^)2Tr,|80/7T, 

7T/180 


. ' 

} 

SELECT DESIRED 
PLANETARY 
CONSTANT UNITS 
FROM TABLE. 




RETURN j 

Figure 6.- CONST flow chart. 
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3.3.5 INPUT (Subroutine) 


Purpose ; Read the input data from the NAMELIST statement, set the default 
values, and initialize all required COMMON block variables. 

Calling sequenoe ; GALL INPUT ($20) 

Called by ; MAIN 

Subroutinea/funotions used ; AEIXYZ, CONST, OUTPUT, XTOPS 

Named COMMON: /CARTC / X( 6 ) ,TIME, ENERGY, R 

/CBASIC / PI,TWOPI,DEG,RAD,DAY 

/CPRINT / PRINT, IPRINT,IELPRT,IUNITS 

/DRAG / CD,AREA,XMASS 

/END / STOP , ISTOP 

/EPOCH / DATE(5),XJDATE 

/KEPLER / EL (6) 

/PERTURB/ IDRAG.ILONG 

/TESS / NMAX,MMAX 

NAMELIST statements: /INPUT/ EL, lEL, STOP, ISTOP, PRINT, IPRTNT, DATE, IDRAG, CD, 

AREA , XMASS , TTESS , NMAX , MMAX , ILONG , lELPRT , lUNITS 

Program data : Size = SOIg f?93io^ words compiled 


FORTRAN 

variable 

Dimension 

Type 

Input/ 

output 

Units 

Description 


AREA 

1 

DP 

I 

2 

m 

Cross sectional surface 
of the satellite ^ valid 
if IDRAG > 1) 

area 

only 

CD 

1 

DP 

I 

— 

Drag coefficient of the satel- 
lite ''valid onlv if TDRAG > 1l 

DATF,^ 

n 

DP 

1 


Cal endar date of epoch 



nATFM) - vf>nr 

(?) - day number 
f3) = hours 
( ))) = minutes 
fb’i r seconds 


f\ 

To tie implemented. 


PS 
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FORTRAN Input/ 

variable Dimension Type output 


Units 


Description 


EL 

fi 

DP 

I/O 

see 

appendix A 

Initial conditions of the satel- 
lite given in Keplerian elements 
or Cartesian coordinates; on 
output it win contain the 
Keplerian elements 






ELm = X or a 

(P) = Y or p. 

(?) = 7 or i 

( 4 ) = X or 'i' 

(■o) = Y or n 

(6) = Y or M 

ENERGY 

1 

DP 

0 

see 

appendix A 

Total energy of the satellite; 
Initially set to n 

TDRAG® 

1 

I 

I 

— 

Flag to determine if drag cal- 
culations are to he included 






= 0 no 
= 1 yes 

I EL 

1 

I 

I 


Flag to determine if the input 
values of EL are given as 
Keplerian elements or Cartesian 
coordinates 






= 1 Keplerian 
= ? Cartesian 

TELPRT 

1 

I 

I 

* 

Flag to determine if the PS 
elements are to be included 
with all output 






= 0 no 
= 1 yes 

TI.ONG 

1 

I 

I 

— 

Flag to determine if the long 
period terms are to be included 






= 0 no fP body orbi t ) 
= 1 no orbit) 

= ? yes^ 

^To be 

implemented . 
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FORTRAN Input/ 

variable Dlmenaion Type output Units Desoription 


IPRINT 

1 

I 

I 


Flag to deter e if the inter- 

mediate printout is to be done 
at a PRINT value of days or 
revolutions 






= 0 no intermediate printout 
= 1 days 
= 2 revolutions 

ISTOP 

1 

I 

I 

— 

Flag to determine if the STOP 
condition is d;./s or revolutions 






- 1 days 
= 2 revolutions 

lUNITS 

1 

I 

I 

— 

Flag to determine what calcula- 
tion constants are to be used 






= 1 kra,seo 
= 2 nm,sec 
= 3 ft, sec 
- 4 m,sec 
= 5 km, hr 
= 6 nm,hr 
= 7 E.r . ,min 

MMAX^ 

1 

I 

I 

— 

Maximum number of tesseral terms 
to be included (valid only if 
ILONG > 1) 

NMAX^ 

1 

I 

I 


Maximum number of zonal terms 
to be included (valid only if 
ILONG > 1) 

PRINT 

1 

DP 

I 

days or 
revs 

Increment at which the interme- 
diate printout is desired (valid 
only if IPRINT > 1) 

RAD 

1 

DP 

I 

rad/deg 

it/ 1 80 

STOP 

1 

DP 

I 

days or 
revs 

Final stop value at which output 
is desired 

TIME 

1 

DP 

0 

hra or 
min or 
sec 

Physical time; initially set 
to 0 

^At 

present 

only MMAX 

= 0 

and NMAX = 2 

are valid inputs. 
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FORTRAN 

variable 

X 

XMASS 


Input/ 

Dimension Type output Units 


Deaoription 


6 DP 0 see Initial state vector 

appendix A 

XCl) ->■ X(3) = X 
X(iO 1. X(6) s V 


1 DP I kg 


Initial mass of the ’.atellite 
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BF^'IN 




Figure 7." INPUT flow chart. 
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3.3*6 OUTPUT (Subroutine’) 


Purpose : Print all desired output during the execution of the ASOP program; 

it contains all output formats and does all unit conversions required 
for output. 

Calling sequence ; CALL OUTPUT (IFORM) 

Called by ; INPUT, MAIN 

Subroutines/functions used ; XYZAEI 

Named COMMON: /CARTC / X(f^l .TIME, ENERGY, R 

/CBASIC/ PI,TWOPI,DEG,RAD,DAY 

/CPRINT/ PRINT, IP, lELPRT.IU 

/END / STOP , ISTOP 

/KEPLER/ XKEP(fi) 

/PS / SIGfUl ,RH0(4),TAU 

Program data : Size = 5738 (379io^ words compiled. 

execution time = ?.8 milliseconds to 33 milliseconds depending 
upon the value of lELPRT. 


FORTRAN 

variable 

Dimension 

Type 

Input/ 

output 

Units 

De.scription 

AEI 

6 

DP 

— 

— 

Character array to accompany 
the Keplerian element output 


AEIdl = 

» A=’ 

(?) = 

' E=' 

f3l = 

' 1=' 

nn = 

' OMEGA = 

fol = 

'NODEr' 

(f^) = 

*Mr ' 


ANG 3 DP — — Character array to accompany 

any angular output 

ANOm r 'DEG' 

(?) r- 'RAD' 
m T ' ' 

CHECK 1 DP — — Energy clieck value 

= fRHO(il) - ENERGY l/RH0(J4l 

DAY 1 DP T sec/dav or Value to convert days into 

mi n/day or hours, minutes, or seconds, 
hr/ day 
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FORTRAN 


variable 

Dimension 

Type 

output 

Units 

De.'icription 

DAYS 

1 

DP 

— 

days 

Print value (DAYS = TIME/DAY) 

DEG 

1 

DP 

I 

deg/rad 

180/ TT 

DST 

7 

DP 

— 

— 

Character array to accompany 


any distance output. 

DST(1) = 'KM' 

(2) = 'NM' 

(3) = ’FT' 

(4) = 'M ' 

(5) = 'KM' 

(6) = 'NM' 

(7) = 'ER' 


ENERGY 

1 

DP 

I 

■■ 

Total energy of the physical 
system. 

HMS 

4 

DP 



Character array to accompany 
any time output. 

HMSCI) = 'DAYS:' 

(2) = ' HRS;' 

(3) = ' MIN:' 

(4) = ' SEC: ' 

lA 

3 

I 



Character array of blanks and 
asterisks (*) 

IA(1) = 

(?) = ' ' 

(3) = 

lELPRT 

1 

I 

I 


Flag to determine if the PS 
elements arc to be printed 

= 0 no 
= 1 yes 

I FORM 

1 

I 

I 


Flag to deterraiiie if the init 


or final condition message is 
to be printed 


=1 initial condition message 
=2 no message (intermediate 
print) 

= 3 final condiMon message 


ORIGINAL 
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FORTRAN Input/ 

variable Dimension Type output 


Units 


ppsorlptlon 


IP 

1 

I 

T 


Flag to determine what print 
condition is being used 

z 1 days 
= ? revolutions 

I STOP 

1 

I 

I 


Flag to determine the final 
stop condition 

= 1 days 
= P revolutions 

lU 

1 

I 

I 

— 

Pointer to the VEL and DST 
arrays 

IXP 

1 

I 



Pointer to the TA a-^ray 

TFORM X 1 -*■ IXP = IP 
IFORM = R + IXP = TSTOP 

REVS 

1 

DP 


revs 

Total number of revolutions 
predicted 

= TAU/Ptt 

RHO 

If 

DP 

I 

see 

PS elements 

RHOm = p RHOn) = P- 

RHO(?) = p^ RHOnn = 

SIG 


DP 

T 

Sf'P 

PS elements 

STOni = a STCnl = rr^ 

SIG(?1 = SlOfin = rr,| 

TAU 

1 

DP 

r 

rad 

Independent variable of the 
PS element's set 

TIME 

1 

DP 

I 

hrs or 
min or 
rpo 

Physical time 

TWOPT 

1 

DP 

T 

rad 

PTf 



FORTRAN 

variable 

Dimension 

Type 

Input/ 

output 

Units 

Description 

VEL 

7 

DP 


_ _ 

Character array to accompany 


any velocity output 


VELO) 

= 'KM/S' 

(2) 

= 'NM/S' 

(3) 

= 'FT/S' 

(4) 

r *M/S ' 

(5) 

= ’KM/H' 

(6) 

= 'NM/H' 

(7) 

= 'ER/M' 


X 6 DP I see Cartesian state vector 

appendix A 

X(1) + XC3) = X 
X(iJ) -f X(6) = V 

XKEP 6 DP I/O see Keplerian elements 

appendix A 

XKEP(l) = a 

(2) = e 

(3) = i 

(4) = u 

(5) = Cl 

(6) = M 

XYZ 6 DP — — Character array to accompany 

the Cartesian state vector X 


XYZ(1) r 

’ X=' 

(2) = 

' Y=' 

(3) = 

' Z=' 

(4) = 

'VX=' 

(5) = 

*VY=’ 

(6) = 

'VZ=' 




OTUC^^-- 

ov 
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Figure 8,- OUTPUT flow chart. 
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3.3.7 POT J2 (Subroutine) 


Purpose ; Compute the potential due to the Jg term of the central body 
geopotential 

Calling sequence ; CALL P0TJ2 (POT) 

Called by ; PSTOX, XTOPS 

Subroutines/funotions us ed; None 

Named COMMON ; /CARTC/ X(8),R 

/KONST/ XMU , XMUI , SQTMU , SQTMUI , EPS 

/RPOOL/ X32,RI,R2I,R3I,XXX(20) 

Program data ; Size = 4 I 5 (33io) words compiled 
execution time = <1 milliseconds 


FORTRAN Input/ 


variable 

Dimension 

Type 

output 

Units 

Description 

EPS 

1 

DP 

I 

see 

appendix A 

= 3/2 y J_ R ^ 

2 e 

POT 

1 

DP 

0 

— 

Potential cue to 

R 

1 

DP 

I 

see 

appendix A 

Magnitude of the position 
vector 

R2I 

1 

DP 

— 

— 

1 /R^ 

R3I 

1 

DP 

— 

— 

1 /R^ 

RI 

1 

DP 

— 

— 

1/R 

X 

8 

DP 

I 

see 

appendix A 

Cartesian state vector 

X(1 ) > X(3) = X 
X(4) ^ X( 6 ) = V 
X(7) = time 
X( 8 ) = total energy 

X32 

1 

DP 




= T? - X(3)^ 
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Figure 9.- P0TJ2 flow chart. 
3t> 




3.3.8 PSANS (Subroutine) 


Purpose : Analytical theory of the ASOP program; it contains only the short 

period and the first order secular terms. 

Calling sequenoe : CALL PfiANS (NN) 

Called by : ASOP, TIMEP3 

Subroutines/functions used ; None 

Named COMMON ; /CBASIC/ PI,TWOPI,DEG,RAD,DAY 

/IPOOL / IMARK 

/KONST / XMU,XMUI,SQTMU,SQTMUI,EPS 

/PS / SIG(11),RH0C4),TAU 

/PSANSV/ FACT0R(4),SIGINI(8) 

/PSTIME/ CLO,FAKTPS,yyy 

Equivalences ; (LS,SIG(4)), (L,RH0(4)), (PHI ,RHO( 1 ) ) , (DSF( 1) ,W( 1 ) ) , 

(DSB(1),W(5)), (S1(1),GAM3(1)), (Y( 1) ,GAM2( 1) ) , (Q( 1) ,GAM( 1 ) ) , 
(HClD.DELBO)), (GC(1),DEL2(D), (P( 1 ) ,DEL( 1 ) ) 

Pr ogram data : Size = 1770s (IOIBiq) words compiled 

execution time = 1 millisecond 

A description of the mathematical symbols used and their relationship to one 
another is given in reference 6; a brief description can also be found in 
Appendixes E and F. Therefore, only a brief mathematical description will be 
given. 


FORTRAN Input/ 


variable 

Dimen-sion 

Type 

output 

Description 

BS 

1 

DP 

— 

1 - H/G 

BY-j 

1 

DP 

— 

1/3 

DY6 

1 

DP 

— 

Ub 

C 

8 

DP 

— 

3*^' j " 1)2, 3, 4 

CN 

1 

DP 

-- 

f1 depending on the value of NN 

COSFC? 

1 

DP 

— 

cos (Ap.j) (see appendix F) 

CnSFC3 

1 

DP 

— 

cc'S 'A^f) (see appendix F) 

CS 

1 

DP 

— 

1/PiG + in = small 'C 

UEL1 

8 

DP 

— 

36/3oj^. 36,''dpj^ k = 1.2, 3, 4 
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FORTRAN 

variable 


Dimenaion Type 


Input/ 

output 



DEL2 

8 

DP 

— 

- 1|2,3t4 

DEL3 

8 

DP 

— 

353 / 30 ^^* 953/3pk ^ = lt2,3,4 

DELTA1 

1 

DP 

— 


DELTA2 

1 

DP 

— 


DELTA3 

1 

DP 

— 

^3 

DIFF2 

1 

DP 

— 

2 2 
P 2 ^2 

DIFF3 

1 

DP 

— 

2 2 
P 3 - a 3 

DSB 

4 

DP 

— 

3 b/aB^^ k = 1 , 2 , 3, 4 

DSF 

4 

DP 

— 

af/36k = 1 , 2 , 3, 4 

EPS 

1 

DP 

I 

e = 3/2 (y Jg r/) 

ETA1 

1 

DP 

— 

sin 0 ^ 

BTA2 

1 

DP 

— 

Hin 20 ^ 

ETA 3 

1 

DP 

— 

sin 30 ^ 

FACTOR 

4 

DP 

— — 

Derivatives of the DS Hamiltonian 
and its combinations (A^, A^, A_ 
and Aj^ in appendix F) 

FAK 

1 

DP 

— 

(2D-3/2 

FAKTPS 

1 

DP 

— 

Derivative of the DS Hamiltonian 
(A^ in appendix F) 

FS 

1 

DP 

— 

f 

FSSQ 

1 

DP 

— 

f^ 

G 

1 

DP 


G 

GAM1 

8 

DP 

— 

3y/3o,^, ay/SPk X = 1 , 2 , 3, 4 

GAM2 

8 

DP 

— 

3Y2^30k> 3Y2/3pji = 1 , 2 , 3, 4 

GAM3 

8 

DP 

— 

3Y3^3o]^> 3Y3'^3pk ^ ~ 1 , 2 , 3, 4 

GAMMA 1 

1 

DP 

— 

1^1 
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FORTRAN 

variable 

Dimension 

Type 

Input/ 

output 

Description 

GAMMA2 

1 

DP 

— 


GAMMAS 

1 

DP 

— 

^3 

GC 

8 

DP 

— 

9G/90j^, 9g/9Pj^ k = 1,2,3,^ 

GIN 

1 

DP 

— 

G-' 

GINSQ 

1 

DP 

— 

g-2 

GM3H 

1 

DP 

— 

G - 3H 

GPH 

1 

DP 

— 

G + H 

GSQ 

1 

DP 

— 


H 

1 

DP 

— 

H 

HC 

8 

DP 

— 

9H/9aj^, 9 h/9Pj^ k = 1 , 2 , 3 ,^ 

HOG 

1 

DP 

— 

H/G 

HSQ 

1 

DP 

— 


IMARK 

1 

I 


Flag determining if one or two 
passes have been made 

0 = 1 pass; 1 = 2 passes 

IQL 

1 

DP 

— 

p/ViT 

L 

1 

DP 

— 

Pj, = ag = L 

LS 

1 

DP 


= A 

NN 

1 

I 


Flag determining if initializing 
or computing 

0 = initializing; 1 = computing 

P 

8 

DP 

— 

9p/9cT^, 9p/3p^^ k = 1,2,3,J| 

P5 

1 

DP 

— 

1/2 

PHI 

T 

DP 

— 

= Og 

P3SQRT 

1 

DP 

— 

yfps 

Q 

8 

DP 

— 

9q/%, 3q/3P^ k = 1,2,3,H 
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PGRTRAN 

variable 

Dimension 

Type 

Input/ 

output 

Description 

QC 

1 

DP 


Q 

QCIN 

1 

DP 

— 

Q"'' 

QCSQ 

1 

DP 

— 


QCV 

8 

DP 

— 

aQ/goj^, 3Q/3p^^ k = 1 . 2 , 3, 4 

QS 

1 

DP 

— 

q 

RHO 


DP 

I/O 

(see SIG) 

S 

8 

DP 

— 

3s/3oj^, 3s/9pj^ k = 1,2,3,11 

31 

i| 

DP 

— 

Derivatives of the generating 
function 

SIG 

4 

DP 


note the location of 

SIG and RHO in COMMON. This 
location makes the equivalence 

= Og, • • • ,pil = ag 

SIGINI 

8 

DP 

0 

The initial values of the a’s 
and p's 

SINPC2 

1 

DP 

— 

sin (Agt) (see appendix F) 

SINFC3 

1 

DP 

— 

sin (A^t) (see appendix F) 

SQTMU 

1 

DP 

I 

yJV 

SQTMUl 

1 

DP 

I 

- 1/2 

y 

S3 

1 

DP 

— 

s 

SUM2 

1 

DP 

— 

1/2 (o 2 ^ + 

SUM 3 

1 

DP 

— 

1/2 ( 03 ^ + P 3 ^) 

TAU 

1 

DP 

I/O 

independent variable j 

TWO 3 

1 

DP 

— 

2/3 

W 

B 

DP 

— 

3w/3aj,i 3w/0pj^ k - 1,2,3,i1 

ws 

1 

DP 

— 

w 

XMU 

1 

DP 

I 

y 





FORTRAN 

variable 

Dimension 

T^Pe 

Input/ 

output 


Description 

XMUI 

1 

DP 

I 

-1 


y 

8 

DP 

— 

9y/9aj^, 

ay/ap,^ k = i 

YS 

1 

DP 

— 

y 


ZET1 

1 

DP 

— 

t-os a-] 


ZET2 

1 

DP 

— 

cos 2a ^ 


ZET3 

1 

DP 

— 

cos 3a ^ 
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BEGIN 





G.HjbjP jq 

and f 





Figure 10,- PSANS flow chart 
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compute: 

Q ,s and c 


I CLEAR 
I WORKING 


array area 


compute: 

8j j Sg ) 

and 


I 

[compute; 


derivatives of 1 

QjSjC j G*-3H*j 


'' 


j COMPUTE: ' 

1 w,y and their, 
derivatives 


\ 

I CLEAR I 

j WORKING 
■ ARRAY AREA j 

1 , 

B 


Figure 10 Continued . 
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3.3.9 PSTOX (Subroutine) 


Purpose ; Transform the^PS (Poincare-Similar) elements into the Cartesian 

coordinates . The subroutines will also oompute the physical 

time for the time iteration stopping procedure. 

Calling sequence : CALL PSTOX (ITIMK) 

Called by ; ASOP 

Subroutines/functions used; P0TJ2 


Named COMMON ; /CARTC / 
/KONST / 
/PS / 
/PSTIME/ 
/RPOOL / 


XI, X2,X3, VI, V2,V3, TIME, ENERGY, R 
XMU , XMUI , SQTMU , SQTMUI , EPS 
srG(i|),RHO(A),TAU 
CLO,FAKPS,Q 

SUMP , T2L , SQT2LI , CS1G 1 , SS1G1 , QCAP , ECOSPH , ESINPH , P , 

ROP , EROOT , EMINPH , FAKT , GC , GCIN , RCAP , SUM3 , X3R0QT , ROOT , 
RCAPDT,XXX(i|) 


Program data ; Size = 3633 (2M3iq) words compiled 
execution time = <1 millisecond 


FORTRAN 

variable 

Dimension 

Type 

Input/ 

output 

Units 

Description 

CS1G1 

1 

DP 

— 

— 

cos 

ECOSPH 

1 

DP 

— 

— 

E cos 

EMINPH 

1 

DP 

— 

— 

E - 

ENERGY 

1 

DP 

0 

— 

total energy of the system 

EROOT 

1 

DP 

— 

— 

p/M 

ESINPH 

1 

DP 

— 

— 

E sin *}> 

FAKT 

1 

DP 

— 

— 

(2D-3/2 

2 ? ^ 

GC 

1 

DP 

-- 

— 

- 1/2 + p^‘ ) = G 

GCIN 

1 

DP 

— 

— 

G"^ 

{TIME 

1 

I 

I 

-- 

Fi.ag determining if only the 
phy.iical time In to be computed 


0 = no; 1 = yes 
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FORTRAN 



Input/ 




variable 

Dimension 

Type 

output 

Units 


Description 

P 

1 

DP 

— 


dist. 


p (see section 2.2) 

P5 

1 

DP 

— 


— 


1/2 

POT 

1 

DP 

I 


— 


Value of the Earth's potential 
due to 

Q 

1 

DP 

0 


— 


q (see section 2.2) 

QCAP 

1 

DP 

— 


— 


Q (see section 2.2) 

R 

1 

DP 

0 


see 


Magnitude of the position 






appendix 

A 

vector 

RCAP 

1 

DP 

— 


— 


R (see section 2.2) 

RCAPDT 

1 

DP 

— 


— 


9R/9t (see section 2.2) 

RDOT 

1 

DP 

— 


— 


Magnitude of the velocity 
vector 

RHO 

4 

DP 

I 


— 


Pl’---’^4 

RI 

1 

DP 



see 


Inverse magnitude of the 






appendix 

A 

position vector 

ROP 

1 

DP 

— 


— 


r/P (see section 2.2) 

SIG 

4 

DP 

I 


— 



SQT2LI 

1 

DP 

— 


— 


P/ V 2L (see section 2.2) 

SSIG1 

1 

DP 

— 


— 


sin 

SUM2 

1 

DP 

— 


— 


1/2 

SUM3 

1 

DP 

— 


— 


0 ^ . p 

3 ^3 

T2L 

1 

DP 



— 



TIME 

1 

DP 

0 


time 


t 

VI 

1 

DP 

0 




Components of the velocity 
vector 

V2 

1 

DP 

0 


app^;nJ 1 X 
A 


VI = X component 

V3 

1 

DP 

u 




V2 = Y component 
V3 ~ 7. component 
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FORTRAN 

variable 

Diir.ension 

Type 

Input/ 

output 

Units 

Description 

XI 

1 

DP 

° i 



Components of the position 




1 


see 

vector 

X2 

1 

DP 

0 


appendix 







A 

XI = X component 

X3 

1 

DP 

0 



X2 = Y component 







X3 = Z component 

X3ROOT 

1 

DP 

— 



- p ^/G 







(see section 2.2) 

XMU 

1 

DP 

I 


see 

P 






appendix A 


XMUI 

I 

DP 

I 


see 

-1 

P 






appendix A 
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COMPUTE- 

ENERGY 


— 

( RETURN ) 

__ 

Figure 11.- PSTOX flow chart. 
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3.3-10 'I’lMEPS (Subroutine) 


Purpose ; Iteration procedure to stop the PS elements at a desired value of 
the physical time. 

Calling sequence ; CALL TIMEPS (TFIN) 

Called by; ASOP 


Subroutines/functions used; PSANS,PSTOX 


Named COMMON: 


/CARTC / X( 6 ),TIME, ENERGY, R 

/CBASIC/ PI,TWOPI,DBG,RAD,DAY 

/KONST / XMU,XMUI,SQTMU,SQTMUI,EPS 

/PS / SIG(8),TAU 

/PSTIME/ CLO.FAKT.QS 


Program data ; 


Size 2026 (130 iq) words compiled 

execution time = 4 to 9 milliseconds depending upon the orbital 
eccentricity. 


FORTRAN 

variable 

Dimension 

Type 

Input/ 

output 

Units 

Description 

CLO 

1 


DP 

I 

time 

Initial value of 

DAY 

1 


DP 

I 

seo/day or 
min/day or 
hr/day 

Value to convert days into 
hours, minutes, or seconds 

DAYS 

1 


DP 

— 

days 

Total number of days elapsed 

FAKT 

1 


DP 

I 

— 

FAKTPS from PSANS = A(, (see 
appendix F) 

lERR 

1 


I 

— 

— 

Error counter 

riEH 

1 


I 

— 

— 

Total number' of irerat. inns 
al 1 owed 

QS 

1 


DP 

r 

— 

q (see section 2.2) 

R 

1 


DP 

T 

see 

appomlix A 

Magnitude of the position vr-i 

TA’J 



DP 

0 

r.id 

Nfrw value of tte independetit 
varlabl'^ {■; ) 

TKIN 

1 


DP 

I 

3-5C or 
min or 
hrs 

Final time desired for .itopp 
t,he i terat ujr. 



1 ,'v 


M9 




FORTRAN 

variable 

Dimension 

Type 

Input/ 

output 

Units 

Description 

TIME 

1 

DP 

I 

sec or 
rain or 
hrs 

Computed value of the physical 
time 

TOL 

1 

DP 

— 

time 

Allowable TOLerance between 
TFIN and TIME that will stop 
the iteration 






|TFIN - TIME] < TOL 

TWO PI 

1 

DP 

I 

rad 

2it 
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3.3.11 XT0P5 (Subroutine) 


Purpose : Transform the Cartesian coordinates (X,V) into the PS (Poincare- 

Similar) elements (u ip) 

Calling sequence : CALL XTOPS 

Called by : ASOP, INPUT 


Subroutines/ functions used: POTJ2 


Named COMMON : /CARTC/ 

/KONST/ 
/PS / 
/RPOOL/ 


X1,X2,X3,V1 ,V2,V3, TIME, ENERGY, R 
XMU , XMUI , SQTMU , SQTMUI , EPS 
SIG(M),RHO(!0,TAU 

G1 SQ , G2SQ , SQT2LI , ESINPH , EROOT , EMINPH , SQTGHI , GCIN , GCSQT , 
SUM3 , R CAP , COS 1 G 1 , SI S 1 G 1 , SUM2 , P , Q , QCAP , ROOT , Z CAP 1 , Z CAP2 , 
ECOSPH,XXX(3) 


Equivalence ; (HC,G3), (PHIC,RHO( 1 ) ) , (LC,RHO(lO) 

Program data ; Size = 5250 (339iq) words compiled 
execution time = <1 milliseconds 


FORTRAN Input/ 

variable Dimension Type output Units Description 

All equations are described in section 2,1 


COSIGl 

1 

DP 

— 

— 

cos 

ECOSPH 

1 

DP 

— 

rad 

E cos 

4> 

EMINPH 

1 

DP 

— 

rad 

E - $ 


EROOT 

1 

DP 

— 

rad 

Vi - 

2Q($ - G) 

ESINPH 

1 

DP 

— 

rad 

E sin 


G1 

1 

DP 

— 

dist2/time 

YV - 
z 

= ‘^x 

G1SQ 

1 

DP 

— 

— 



G2 

1 

DP 

— 

dist^/time 

ZV - 

X 

""z = 

G2SQ 

1 

DP 

— 

— 

0 

G 

y 


G3 

1 

DP 

— 

dist'^/time 

XV - 

V 

= G, 
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FORTRAN Input/ 


variable 

Dimension 

Type 

output 

Units 

Description 

GC 

1 

DP 

— 

dist^/time 

g = ,/g^ + g^ + g^ 
V X y z 

GCIN 

1 

$ 

DP 

— 

time/dist^ 

1/G 

GCSQ 

1 

DP 

— 

— 

0^ 

GCSQT 

1 

DP 

— 

— 

VF 

HC 

1 

DP 

— 

— 

H = =3 

LC 

1 

DP 

0 

(dist/time)^ 

L = = Og 

P 

1 

DP 

— 

dist 

p (small p) 

PHIC 

1 

DP 

0 

dist^/tirae 

5> = P, = Pg 

POT 

1 

DP 

I 

— 

Potential due to 

Q 

1 

DP 

“ 

— 

(i - “) 

QCAP 

1 

DP 

— 

— 

Q 

R 

1 

DP 

0 

see 

appendix A 

The magnitude of the position 
vector (r) 

RCAP 

1 

DP 

— 

— 

R 

ROOT 

1 

DP 

— 

— 

3R/at 

RHO 

i} 

DP 

0 

— 

^ in RHOd) ♦ RHOCd 

RT 

1 

DP 

— 


Inverse magnitude of the posi' 
vector magnitude ( 1/r) 

RSQ 

1 

DP 

— 

-- 

Square of the position vector 
(r") 

SIG 

i) 

DP 

0 

see 
ref. 1 

-*• in SIG(l) SIG(4) 

SIS1G1 

1 

DP 

— 


sin 

GQT2LI 

1 

DP 

— 

— 

u/ 

GQTGHI 

1 

DP 

— 

— 

-\J2/\.G + H) 
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FORTRAN Input/ 

variable Dimension Type output Units Desoription 


SUM2 

1 

DP 

— 

dist2/tirae 

$ _ G 

SUM3 

1 

DP 

— 

dist^/tirae 

2{G - H) 

TAU 

1 

DP 

0 

rad 

Independent variable; initially 
set to zero (0) 

TIME 

1 

DP 

I 

time 

Initial physical time 

VI 

1 

DP 

0 

see 

Components of the velocity 
vector 

V2 

1 

DP 

0 

‘ appendix 

A 

VI = X component 

V3 

1 

DP 

0 


V2 = Y component 
V3 = Z component 

XI 

1 

DP 

0 

see 

Components of the position 
vector 

X2 

1 

DP 

0 

^ appendix 

A 

XI = X component 

X3 

1 

DP 

0 


X2 = Y component 
X3 = Z component 

XMU 

1 

DP 

I ' 

see 

appendix 

y 

XMUI 

1 

DP 

I . 

A 

1/U 

ZCAP1 

1 

DP 

— 

— 

^1 

ZCAP2 

1 

DP 

— 

— 

Zg 


5^1 



compute: 

ANGULAR 
MOMENTUM 
VECTOR S' 


i 

(Fall potj2 

COMPUTE: 

La 



I 


CALCULATIONS 

FOR 


CALCULATIONS 

FOR 



CALCULATIONS 

FOR 



Figure 13.-XT0PS flow chart. 
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3.3.12 XYZAEI (Subroutine) 


Purpose i Transform the Cartesian coordinates (X,V) into the Keplerian elements 
(a,e,i,o),S^,M) 

Calling sequence ; CALL XYZAPI 
Called by ; OUTPUT 
Subroutines/ functions used: None 


Named COMMON ; /CARTC / 
/CBASIC/ 
/KEPLER/ 
/KONST / 
/RPOOL / 


XI , X2,X3, VI, V2,V3, TIME, ENERGY, R 
PI,TWOPI,DEG,RAD,DAY 
A,E,XI,OMEGA,XNODE,XM 
XMU , XMUI , SQTMU , SQTMUI , EPS 

VSQ , RRDOT , ECOSE , E3INE , P , H , WX , WY , WZ , EA ,RCOSF , RSINF , 
COSWF , SINWF , TEMP , RCOSL , RSINL , XXX ( 7 ) 


Program data : Size 4068 (262 iq) words compiled 

execution time = <1 railliseoonds 
Considers only elliptic motion 


FORTRAN 

variable 

Dimension 

Type 

Input/ 

output 

Units 


Description 

A 

1 

DP 

0 

see 

appendix 

A 

Seraimajor axis of the orbit 

COSWF 

1 

DP 

— 

— 


Cos ('jJ + f) 







A) = argument of perigee 
f = true anomaly 

E 

1 

DP 

0 

— 


Orbit eccentricity (e) 

EA 

1 

DP 

— 

rad 


Eccentric anomaly (E) 

ECOSE 

1 

DP 

— 

rad 


E cos e 

ESINE 

1 

DP 

— 

rad 


E sin e 

H 

1 

DP 

— 

dj st^/tirae 

Total angular momentum 

OMEGA 

1 

DP 

0 

rad 


Argument of pericenter (w) 

P 

1 

DP 

— 

-- 


1 - e 

R 

1 

DP 

0 

see 

appendix 

A 

Magnitude of the position 
vector 
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FORTRAN 



Input/ 



variable 

Dimension 

Type 

output 

Units 

DescripL i on 

RCOSF 

1 

DP 

— 


name as R 

r cos f 







f : true anomaly 
r : R 

hCOSL 

1 

DP 

— 


same as R 

r cos L 







L = mean anomaly 
r = R 

RRDOT 

1 

DP 




see 

X • V 






appendix A 


RSINF 

1 

DP 

— 


same as H 

r si n f 







r = R 

f = true anomaly 

RSINL 

1 

DP 

— 


same as R 

r sin L 







r = R 

[. = mean anomaly 

SINWF 

1 

DP 

— 


— 

sin (to + f) 







iu = argument of perigee 
f = true anomaly 

VI 

1 

DP 

I 


see 

Components of the velocity 
vector 

V2 

1 

DP 

I 

) 

append! x 
A 

VI = X component 

V3 

1 

DP 

! 



V2 = y component 
V3 = 2 component 

VSQ 

1 

DP 

— 


see 

Magnitude of the velocity 
veclor, squared, (V^) 






appendix A 

wx 

1 

DP 

— 


di st^/'time 

Components of the total 
anguiar momentum 

WY 

1 

DP 

-- 

' 

WX = X component 


WZ 

1 

DP 




WY = Y component 
W7. = Z component 
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FOKVHAN Input/ 

variable Dimension Type output 


Units 


Descr ipclon 


XT 

1 

DP 

I 

see 

Components of the position 
vector 

X2 

1 

DP 

I 

appendix 

A 

X1 = X component 

X3 

1 

DP 

I 


X2 = y component 
X3 = Z component 

XM 

1 

DP 

0 

rad 

Mean anomaly 

XMU 

1 

DP 

I 

see 

appendix A 

Central body gravi tatiorial 
constant (y). 

XHODE 

1 

DP 

0 

rad 

Argument 0 )' the ascending 
node (SI) 
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Page 1 of 2 


Figure 14.- XYZAEI flow chart. 
59 


05' POOK quality 










Figure T4.- Concluded, 


Page 2 of 2 
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3.^1 LABELED COMMON; 


Notation: R - Real variable | 

I - Integer variable i Type 



S - Single 
D - Double 

precision 1 
precision | 

Precision 


/CARTC/ 

Looation 

Name 

Dimension 

Type/ 

precision 

Description 

1-3 

X 

3 

R/D 

Cartesian coordinates for the 
position of the satellite; 
X,Y,Z 

lt-6 

V 

3 

R/D 

Velocity vector of the satel- 
lite; V^,Vy,V^ 

7 

TIME 

1 

R/D 

Elapsed time (hr, rain, or 
sec) 

8 

ENERGY 

1 

R/D 

Total orbital energy 

9 

R 

1 

R/D 

Magnitude of the position 
vector 


In subroutines ; MAIN, AEIXYZ, ASOP, INPUT, OUTPUT, 
P0TJ2, PSTOX, TIMEPS, XTOPS, XYZAEI 


/CBASIC/ 



Looation 

Name 

Dimension 

1 

PI 

1 

2 

TWOPI 

1 

3 

DEG 

1 

i) 

RAD 

1 

5 

DAY 

1 



R/D 2n 

R/n tt/180 

R/D l80/?[ 

R/D Converts Jays into hours, min- 
utes, or seconds 


In subroutines ; MAIN, ASOP, CONST, INPUT, OUTPUT, PSANS, 
TIMEPS, XYZAEI 



6 1 



/CPRINT/ 


Type/ 


Location 

Name 

Dimension 

precision 


Description 


1 

PRINT 

1 

R/D 





2 

IPRINT 

1 

I 


Input 

4.1.1 

parameters, see 

sectlc. 

3 

lELPRT 

1 

I 





ii 

lUNITS 

1 

I 





In subroutines: MAIN, 

CONST, INPUT, 

OUTPUT 





/DRAG/ 



Type/ 





Location 

Name 

Dimension 

precision 


Description 


1 

CD 

1 

B/D ’ 


Input 

parameters, see 

section 

2 

AREA 

1 

R/D 


4.1.1 



3 

XMASS 

1 

R/D 





In subroutines: INPUl 







/n:ND/ 



Type/ 





Location 

Name 

Dimension 

precision 


Description 


1 

STOP 

1 

R/D 


Input 

4.1.1 

parameters, see 

section 

2 

ISTOP 

1 

I 





In subroutines; MAIN, 

INPUT, OUTPUT 






/EPOCH/ 



Type/ 





Location 

Name 

Dimension 

precision 


Description 


1-5 

DATE 

5 

R/l) 


Input 

4.1,1 

parameter, nee 

section 

6 

XJDATE 

1 

R/D 


Julian date 



In subroutines; INPUT 
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/GEO/ 


Location 

Name 

Dimension 

Type/ 

precision 

Description 

1 

RE 

1 

R/D 

Equatorial radius of the cen- 
tral body 

2 

CJ2 

1 

R/D 

Jp geopotential coefficient 
or the central body 

In subroutines: CONST 




/IPOOL/ 

Location 

Name 

Dimension 

Type/ 

precision 

Description 

1-10 

— 

— 

I 

Temporary integer variables; 
this COMMON block is used to 
help save data storage within 
the ASOP program. 

In subrou*' 

ines: PSANS 




/KEPLER/ 

Location 

Name 

Dimension 

Type/ 

precision 

Descriptinn 

1 

A 

1 

R/D 

Semiraajor axis of the orbit 
(a) 

2 

E 

1 

R/D 

Eccentricity (el 

3 

XI 

1 

R/D 

Orbital inclination to the 
Equator ( i ) 

i| 

OMEUA 

1 

R/D 

Argument of pericenter (..)) 


XNODE 

1 

R/U 

Argument of the ascending (it) 
node 

6 

XM 

1 

B/D 

Mean anomaly (M) 

[n subroutines: AEIXYZ: 

, INPUT, OUTPUT, XYZAEI 





f)3 

ORIGINAL PAGE IS 
OF POOR QUALTIY, 



/KONST/ 


Type/ 


Location 

Name 

Dimension 

precision 

Description 

1 

XMU 

1 

R/D 

Gravitational constant of 
the central body (M/ 

2 

XMUI 

1 

R/D 

1/y 

3 

SQTMU 

1 

R/D 


4 

SQTMUI 

1 

R/D 

i/Vu 

5 

EPS 

1 

R/D 

2 

£ = 3/2 M J., R^ where 

y ft 


Jg = J 2 geopotentlai 
coefficient 


R = Equatorial radius of 
the central bndy 

In subroutines : AEIXYZ, CONST, P0TJ2, PSANS, 

PSTOX, TIMEPS, XTOPS, XYZAEI 


/PERTRB/ 



Type/ 



Location 

Name 

Dimension 

pj eels ion 

Description 


1 

IDRAG 

1 

I ) 

Input parameters, see sec- 




\ 

tion 4.1.1 


2 

ILONG 

1 

I ) 



In subroutines: INPUT 





/PS/ 



Type/ 



Location 

Name 

Dimension 

precision 

Description 


1-4 

SIG 

4 

R/D 

PS elements ^ 

‘^ 3 , ^4 

5-8 

RHO 

4 

R/D 

PS elements P^, P^, 

, P 4 

9 

TAU 

1 

R/D 

Independent variable 

of the 


PS elements. 

In subroutines : MAIN, ASOP, OUTPUT, PSANS, PSTOX, 

TIMEPS, XTOPS 



/PSTlMhV 


Location 

Name 

Dimension 

Type- 

precision 

Description 

1 

CLC 

1 

R/D 

Initial value of or (set 

when initializing, wEWX = 0) 

2 

FAKTPS 

1 

R/D 

9a ^^/3 t 

3 

Q 

1 

R/D 

q 


In Bubpoubines ; ASOP, PSANS, PSTOX, TIMEPS 


/RPOOL/ 

Location Name Dimension 


Type/ 

precision 


Description 


1 - 2^1 


24 R/D Temporary real variables; 

this COMMON block is used 
to help save storage within 
the ASOP program 


In subroutines ; AEIXYZ, P0TJ2, PSANS, PSTOX, 
XT0P3, XYZAEI 


/TESS/ 

Location 


Type/ 

Name Dimension precision 


Descri < "in 


1 NMAX 

2 MMAX 


1 

1 




Input pararaetf, , see sec- 
tion ^.1.1 


In subroutines: CONST, INPUT 
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4.0 USER *3 GUIDE 

This aeotlon is intendsd to give the user all the information neoessary to oper- 
ate the ASOP programs. Beoauae the program is designed to operate in two modes 
{stand-alone and subroutine package), eaoh mode of operation is described 
separately. 

The first part of this section (section 4.1) will describe the general input 
parameters and options available when using the stand-alone ASOP program. Also 
described in this section are the standard default values and the typical con- 
trol cards needed for executing the program lu the demand or the batch mode. 
Finally some sample output is given to help the user if modifications are to 
be made to the system. 

Section 4.2 will deal with the ASOP subroutine package. This section will de- 
scribe the necessary modules that are used within the package, as well as any 
interface requirement that the user must be aware of if he is to include this 
package in his own software. The input to and output from the \S0P subrou- 
tine are also fully described in this section, as are the subroutines default 
values . 


4.1 INDEPENDENT PROGRAM 

The ASOP program was designed as an interactive program capable of giving the 
user fast, accurate answers to Shuttle-type orbit problems. The program may 
also be run in a batch environment if a large number of cases must be 
investigated . 

There are two basic methods used to control the operation of the ASOP program; 
flags and direct-user interaction. The flags are used to indicate the type 
of data being entered and to select certain options within the program. Direct 
interaction allows the user to check his data to insure accuracy before contin- 
uing or to terminate the program. 

Primary data input to the ASOP program is accomplished using the NAMELIST 
*$INPUT'. The necessary input variables and user responses to program questions 
are described in section 4.1.1, and the program default values are described 
in section 4,1.2. Section 4.1.3 explains the printed output generated by the 
ASOP program. Finally, sections 4.1.4 and 4.1.5 describe the instructions ne- 
cessary to run the ASOP program and give an example of the resulting output. 


4.1,1 Input Description 

The NAMELIST is the primary method of getting data into the ASOP program. How- 
ever, during normal operation, the user lo expected to interact with the program 
by supplying additional information. After starting the ASOP program (see sec- 
tion 4.1.4), it will ask for the NAMELIST data with the statement 


INPUT DATA USING NAMELIST '$INPUT* 
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At this point, the user has three options: 

a, To enter the NAMELIST data directly from the keyboard. 

b, To add a data file containing the NAMELI.ST information using the SADD com- 
mand (ref, 5). 

c, To enter @E0F to terminate program execution. 

If option C is selected, the program will respond with 

**N0RMAL PROGRAM TERMINATION** 


and the program will stop. If option A or B is selected, than the program will 
print out all the NAMELIST variables and their associated values (including de- 
fault), as well as the initial conditions of the problem. The message 

ENTER: C = CONTINUE; S = STOP 


should then appear. Here, the user should check the input data and enter the 
necessary letter. If a C is entered, the program will continue the execu- 
tion as directed by the input. When the input stop condition is satisfied, the 
program will again ask for data input as described earlier. The series of in- 
structions can be repeated as often as necessary. 

If an 3 is entered in response to the message, then the program will respond 
with **NORMAL PROGRAM TERMINATION**. The user may then make the necessary cor- 
rections to the data and restart the program as described in section 

Table X describes the input variables that may be used in the NAMELIST '$INPUT'. 
Whether keying in the information or creating a data element, a ^$INPUT must 
be entered first where the X represents one or irore space.s. Each variable 
entered must be preceded by one or more spaces, and if more than one variable 
Is to appear on a line, they must be separated by a blank or a comma (,). To 
terminate the NAMELIST input a )i($END or must be the last item entered. 

(See ref. 7, pages fi through 13, for a complete description of a NAMELIST 
statement) . 
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TABLE I.- NAMELIST INPUT VARIABLES 


Variable 

Type 

No. of inputs 
required 

EL 

DP 

6 


lEL I 1 


STOP DP 1 


ISTOP I 1 


PRINT DP 1 


& 

Default value. 


Description and 

available options 

Must be supplied by the user. Can be 
given as Cartesian coordinates or 
Keplerian elements as determined by 


the 

flag 

ILL 


EL 

CO 

— 

X 

or 

a 

EL 

(2) 

— 

Y 

or 

e 

EL 

(3) 

= 

Z 

or 

i 

EL 

(4) 

= 

X 

or 

u 

EL 

(5) 

= 

Y 

or 


EL 

(6) 

— 

Z 

or 

M 


All angular input is given in degrees. 
All other values must correspond to 
the option specified by the flag 
lUNITS. 

Flag determining the type of initial 
conditions input. 

1 = Keplerian elements^ 

2 = Cartesian coordinates 

Final condition that must be satisfied 
in order to stop program execution 
normally 

Flag that specifies the type of STOP 
condition 

1 = STOP in days 

2 = STOP in revolutions^ 

Increment for the printed output. 

A value Is not needed if IPRINT is 
set to 0. PRINT = 0.0 is a valid 
entry. 
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TABLE I,- Continued 


Variable Type 
IPRINT X 


DATE^ 


DP 


IDRAG^ I 


area'^ dp 


XMASS'^ DP 


ILONG I 


No. of inputs Description and 

required available options 

1 Flag that specifies the PRINT increment. 

0 = No PRINT increment^ 

1 = PRINT is days 

2 3 PRINT is revolutions 

6 Date of ep )ch given as a calendar 

date of the form 

Day, month, year, hours, minutes, 
seconds 

(to be implemented) 

1 Flag that specifies if the drag equa- 

tions are to be included in the 
computation 

0 = No^ 

1 = Yes 

(to be implemented) 

1 Frontal surface area of the satellite. 

A value is not needed if IDRAG is 
set to 0 

(to be implemented) 

1 Total mass of the satellite. A value 

is not needed if IDRAG is set to 0 
(to be implemented) 

1 Flag that specifies the type of poten- 

tial terir,.:^ to be included in the 
computations 

0 = None (two-body orbit) 

1 = short period, and first 

2 * D 

order secular terras 

2 = long period, and higher 

order secular terms'^ 


^Default value. 

To be implemented. Current values set to zero. 


ORIGINAL PACE IS 
OF POOR QU/' 
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TABLE I.- Concluded 


Variable Type 
NMAX I 

MMAX I 

lELPRT I 


I UN ITS I 


No. of inputs Description and 

required available options 

1 Total number of zonal terms to be in- 

cluded by the geopobential model; a 
value is needed only if ILONG is 
set to 2 (to be implemented). 

1 Total number of tesseral terms to be 

included by the geopotential model; 
a value is needed only if ILONG is 
set to 2 (to be implemented). 

1 Flag to determine if the PS elements 

are to be included with all printout. 

0 = No^ 

1 = Yes 

1 Flag that specific ' he units of the 

input data and se the appropri- 

ate physical const 

1 = km, sec^ 

2 = nm, sec 

3 = ft, sec 

4 = m, sec 

5 = km, hr 

6 = nm, hr 

7 = E.r., min 


^Default value. 
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4.1.2 Default Values 


To help shorten the number of data values that must be supplied by the user, 
the ASOP program assumes certain default values for those variables not expli 
oitly mentioned on the input NAMELIST. These default values are listed in ta 
ble II and a desoription of the variables can be found in section 4,1.1, Any 
variable not listed in table II must be specified by the user. 


TABLE II.- DEFAULT NAMELIST VALUES 


Variable 

Default value 

I EL 

1 

.rop 

100.0 

I3T0P 

2 

PRINT 

0.0 

IPRINT 

0 

IDRAG^ 

0 

AREA^ 

0.0 

CD^ 

0,0 

XMASS^ 

0.0 

ILONG 

1 

NMAX^ 

2 

MMAX^ 

0 

lELPRT 

0 

lUNITS 

1 


^To be Implemented. 
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^.1.3 Output Desorlptiog 


After the ASOP program has been started with the command 
eXQT »NUMEG. ASOP -PROG 


and the input data has been added, the program will print out all of the 
NAMELIST variables as shown in figure 15 and the initial conditions shown 
in figure 16. 


i 


f INPUT 
EL 

_ 

. 327265861O000G00OGD+0O4 

9SD+00H, 


. OO0OO0O000OOOO0OOOD+00O 

0OD+OG1 , 


. 39 1 02S30 1 999999999D+G0 ] 

OOD+000 



lEL 

- 

+ 2 

CTCiP 


. 1OGGOOOOGOO0OGO00OD + OO3 

I STOP 

i; 

+ •■> 

PRINT 

- 

. 5OOOGOnoO00OOGeOGaD+GG2 

IPRINT 


+ 2 

OPITE 

- 

. C J0GOGeOGGO0O0G0GOD + O0O 

OOD+000, 


. 030GP.GOOGOGOGOGGD-I-OOO 

O0D+0OG, 


. 0OGGGGGOOOGOGO0GOGD-H0OG 

I DRflij 

“ 

+ 0 

RE W 


. GGOOG0GOOOGGGQGOG0D+O00 

CD 


„ O0O0GOGOOOeG0GGOGOD+GOO 

KMflSS 


. 00GOGOGG0O0OOOOOOOD-I-OOO 

ILGNG 


+ 1 

NMRX 


-1-2 

MMRX 


-i-G 

lELPRT 


-1-1 

I UN I TS 
TEND 


+ 1 


. 57-37872519999999 
S77H77591QQO0GOG 
. 0GOGOOOGOOf)0OO0O 


. 0000000000000000 

. 0000000000000000 


Figure 15.- NAMELIST data output format. 


□nv'-;; 0 .. onnoooo 


INM [i'll r i.NniTH.iNS 

! PF VO : 0 . 0000000 ' CHECK : - 1 , OOOOOOO hOO 


R:-.: 6.b993.359+03 KM ' K - l-hor, 
iiMEF.M.: .3.39. 5?fi 'H [iPra iHi'tlTr. .. r ... 

S 272<:..bCt. v0.7 |: M ! T- 5.73 .C 
v::_-6 ,. 77'i ,''?’09 1-00 i h/-- 1 vv:- 3. our: 
[•'S LLf Mr- MT'-;.:: 

1 o'' 2 Hor.:.’t I 001 -. 197.300.34 ^ orn 
' . 1 1. 7', :,:L9i 00^ , 24F7 ’ K 'tP < On ) 


+ 


n . 'i 0 ? ! 

U: . GGOOO 

DEL 

'id hf-e I 

M-:: 19.P3''>7..3 

rc n 

. r. 10 - ': t M I 

Z-' Q.OGOOdOO 

!- M 

n .G0 I fVS! 

VZ= 0. 0O0OGOG 

1 M / 

. ndOOGOiio 

. 89 3 39(.,3f'. + 00 1 


f 'fif 1 d 1 . < 1 0 

. >9 77 /b'-'i 1-00,? 



ER I I (ONriMIlE, C 


Figure 16.- Initial condition format. 


After the initial conditions have been displayed, the program will wait for the 
user to check the input parameters. Some checks that can be made are 
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a. CHECK value should be -1.0000000+00 


b. DAYS and REV3 should be 0.0000000. 

Q. A double asterisk (*'*') will appear ai'cc-* the output condition that is being 
satisfied, i.e., after the DAYS or the REV'S value. This will agree with 
input value of IPRINT or ISTOP. 

d. Are the initial oonditions (a, e, 1, eto.) the desired values? 

e. Are the units oorrect (kra, kw/a , eto.)? 

Onoe the user is satisfied that the printed Initial oonditions are the ones 
needed , a C should be entered . When the C is entered , the program will 
continue execution and print out information as speoified by the input 
parameters, i.e., STOP, ISTOP, PRINT AND IPRINT (figure 17). 


ENTER: 

>C 

*• c • 


C = CONTINUE; S = 5T0P 


! REVS: E. 0000000+01 


+ ! CHECK: 4 . 99OO3S0 ~0: 


A — 

6. 6993637 + 03 

KM 1 E = 

1 . 501 3305-02 


1 = 

. 00000 

DEG 

omedh= 

359 53545 

DEG INODE- 

. 00000 

DEG i 

M- 

. 450 ■>. 

DEG 


2. 3173764+03 

KM ! Y = 

5.9670 

■+77 + 03 

KM 

2- 0. 

OOOOOGO 

KM 

vx= 

-7. 0:302476 + 00 

KM/S! VY= 

3 , 54 3 '7 

'707 + 00 

KM/S 

IV 2 --:; 0 

0000000 

LM, 


PS ELEMENTS: 

. S 1 5 2 S 9 9 '’I' + 0 0 3 
. B1C7HS73+OOE 


. 2fJbE?7H2+Q01 
. 13£;3fc.ES:5 + 0n I 


. 00000000 
. 00000000 


. 2 721 09 2H ^•OQb 
. 29 7 7''': S'-:' + 007 


FlNfU- CONDI T IONS 

OftVS: 6,2937974+00 I REVS: I . OOOOOOO+O; 

6, 6'S933H4 + 0.3 I M I E= J . 49 1 . ;6--02 

0MEEri= 33.64542 fJEG ! N0DE= .00000 DEG 

:<= 2. 35G3<^ ■•1+03 KM ! Y= 6 . 1 7?4332 + 03 KM 
VX=-7. 3242597+00 KM/SI VV~ 2,7434611+00 LIVS 
PS eu:ments: 

.02952542 + 003 3 1 747392 <-00 1 . OOOGGOOG 

. 51b7H349+v.05 -.95472 769-003 . 00000000 

INPUT Dl-irn Li'SlNG NFIMELIST ' 1 INPUT' 


+ + ! CHEC 

I- 
M = 

Z = 

V Z 


I;: 4 . 968'+0:'0 -07 

. OOGOQ DEG 

34 1-060':':': PI.LG 

0.0000000 Mi 

O.OOOOf‘00 I M.'' 


. 54420755+00' 
2 '+; 7 ':' '■ ‘1 + 00 :; 


Figure 17.- Intermediate and final ou'-^put format. 


In general, all output is clearly labeled, but some terras should be explained 
further. 

DAYS; Total number of days elapsed since the starting epoch. 

REVS; Total number of revolutions completed, 

CHECK; Value indicating the accuracy of the analytical theory; although 
this value is necessary as a check on the theory, it la not a 
sufficient check. 


OF POOK QUAUn 
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PS ELEMENTS: The Poinoare-Siinilar elements listed as 



Double asterisk (**): Indicates the stopping condition being satisfied} 

this flag will move between the DAYS and REVS 
value as needed. 

4 . 1 . Run Setup (Control Cards) 

The ASOP program is written in standard FORTRAN V and designed to run on the 
NASA/JSC UNIVAC 1110 computer using the EXEC-8 operating system. All the re- 
locatable and executable elements are on the file FM6-N08569*NUMEQ. ASOP may 
be executed by entering the following for demand operation. 

a. gQUAL ?I'!'6~N08569 

b. @ASG,A *NUMEG. 

n. @XQT »NUMEG.AS0P-PR0G 

d. Add input data 

e. Enter the letter C or S (see section 4.1.3) 

f. §E0F or go to step d. 

g. eFIN 

If run in a batch mode, the following input cards are needed. 

a. §QUAL FM6-N08569 

b. §ASG,A «NUMEG. 

o, @XQT »NUMEG.AS0P-PR0G 

d. Add data file or data cards 

e. 'C' 

f. Repeat instructions d and e as often as peoessary. 

g. @E0F 

h. §FIN 
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^*1.5 aatnpXe Coapufcer Rur) 


In this seoticuit a sample ooraputer run is reproduoed for a typioal shuttle- 
type orbit. The orbit hes been predicted for 100 revoiutians (?r6.3 days) with 
the output given every three days, 

This example jjs intended bo familiarize the user with the format of the ASOP 
output and bo illustrate ';be use of the various input options dieousaed in seo- 
tion A full description of the output is given in .-.eotion ^1.1.3* 


Initial pararoebera 

Serairaajor axis (a) 
Eqoentricity (e) 
Inclination (1) 

Argument of perigee (w) 
Argument of the ascending 
Mean anomaly (M) 


6699.3532 km (1 ,05 BR) 

,001 

30 degrees 
18 degrees 
(fl) 20 degrees 
22 degrees 


Input parameters under the NAMELIST $INP0T 


$INPUT 



EL(1) 

5 

3542.07055 

EL(2) 


5256.17858 

EL(3) 

= 

2152.20751 

EL(4) 

- 

-6.41478310 

EL(6) 

= 

3.11545050 

EL(6) 


2.95692882 

lEL 

= 

2 

STOP 

r 

100,0 

ISTOP 


2 

PRINT 


3.0 

IPRINT 

“ 

1 

ILONG 


1 

NMAX 

= 

2 

MM AX 

"S 

0 

lELPRT 


1 

lONITS 

$END 

X 

1 
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SampTe computer run, 

JfXQT #WUMfft,flSOP-PRO® 

INPUT OflTft USING NRMELIST ’$INPUT’ 
>tfiOO S.0«T0-fiS0P3 

> STOP = IO0.0DO 

> ISTOP = 2 

> PRINT = 3.0DO 

> IPrlNT = 1 

> »EKO 
tINPUT 


EL 

:r 

. 3542070543999999990+004 

99D+00H, 


.2152207510000000000+004 

99D+051 , 


. 3115450500000000000+001 

000+001 



lEL 

“ 

+ 2 

STOP 

z 

. IOOPOOOOOO00OOOOOCD+Or3 

ISTOP 

zz 

+ 2 

PRINT 


. 3OPOGOOOOOOOOOOGOOD+G01 

IPRINT 


+ 1 

DOTE 

OOD+000, 

= 

. OOQOOOeGOOOOOQOOOOD+OQO 

000+000, 


. OO0GOOPOOOOGOOOOOOD+O0O 


* 

IDRfiG 


. OOOOOOOOOOOOOOOOOOD+OOO 


■4* 

RRER 


. GOOO0OOOG00OO00O0OD+0OO 

CO 

— 

. 0400009000000000000+000 

XMRSS 


. oooooooooGoeenooooo+OGo 

ILONG 


+ 1 



+ 2 

MMPK 

— 

+0 

lELPRT 

- 

+ 1 

lUNITS 

#END 


+ 1 


.52E6178S79'3S9939 
6Hm78309'3<^'39999 
. 29r.f:. 929820000000 


.OOOOOO0OOOOCCOOO 

. ooooeoooooooecoo 


INITIRL CONDITIONS 

DRYS: O.COOOOOO i ’J' ! REVS: O.OGOOOOO ! CHEC - 1 . OOOOOOO+OO 

0= S. 6 ;>'V3r;32+-03 KM i £= 9 . 099r.t'.0r.-0H f l~ 70 0O9OQ DFG 

'■■M'rnf!- 17.9;?T6 OfG 'NODE- 20.0000© HtTi ! Hr 22.0S&98 DEH 

:r 3.E.'»19.0*;+0.3 KM * Y= 5.2S61 /::.r.*0 3 f .1 » 7- 2. 1S220fG+O3 KH 

Vr's-C . 'tlH /0.--1+00 KM/S> VY= 3. 1 JiSHEor.+n© VZ* 2, 9'-.6915.0+n0 \ 

rs 1-:L EMFHVS: 

. 10’. 7SH>: fOCl . 6379049N-00? >. 002 . B9.3 KJPT'i ♦■000 

. KU.7SK! '1+C"' .KAIG97’:' -o;.l . ; 1 0<-,?H'i7 Ki,i7 . ;-97»s9'! : 9+002 

I RIT,!.. r r ( nN'miMF , ■?: - STOP 
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Sample computer run,- Concluded 


ENTER: C = CONTINUE; S = STOP 

>C 

nOYS- 3.0000000+00 ++ I REVS: 4.7532597+01 'CHECK: 1.0G331O7-C6 

^ . + 1 

fl:- 6. 596 Hie 01+03 KM ! E= 9.60^^1310-04 ' I:: 29. 9 , C53 KCu 

OMEGH-.: 296.00154 0+G 'NODE; 35t-, jnssO DEG ' .33F.,7:;;775 OF ’, 

Y' = -3. 47.53360+01 KM ! Yi-G. V-)5? 17.3 +0-3 KM » Z--C-. 3422.. TH03 KM 

7.7207336 + 00 1 M/S ! V 57r.fl9.5 1 - 02 KM/3'Vr,- 9 . .'OtmC- 02 1 !• •'■? 

F-. ELEMENTS: 

. .30001249+003 . 1432S777 ' 000 . 3336.''65i l OOl . I .'f' 3 1 tOC6 

.51664463 + 035 , 396*Y 573-00 I ,, 1 U 35 + 00.5 . 2'-:« ; .94 23 i 00 ? 


DRY'S: 6., 0000000 + 00 

* 

\ f VZ 

: 9.5165702+01 

•CHECK: 5. .35 73rr • 

*i ; 

A= 6.6'^9S87C+C3 

KM 

! Ex 

7. 71 '+4 707 -O': 

! Ix 30.00177 

CEl. 

OMEDAx 108-73237 

DEG 

INODE:: 

536, 21553 

DEG I h: 3't.064 19 

PEG 

>:=-3 4707616+03 

KH 

I V- 

5.35655)0+03 

KM I 2:- 2. 021 +'. ,8+03 

y i" 

v';:x-6. 4158 754+00 

KM/C 

' VYx 

“2. 992L'5Cr + 00 

KM/S' V2x-3. 07't9‘. 1.7+00 

KK/S 


V . . “-I? , T 1 / 5T + 00 V f jz- ^ i-. n^':) ’ v i: - ~ . n itv > 

PS ELEMENi'L-.-. 

.59898062+003 -.3507305.3-001 .47459491+002 .51340075+006 

.51676525+005 .94153688-001 .10768413+005 . 29769423+0C2 


FINAL CONDITIONS 

30YS: 6-3047894+00 ! REVS: 1.0000000+02 +*!CHECK: 1.0044341-C6 

+ +-- 

Q- 6.6964889+03 KM ! £= S.4V15012-Q4 ' 1= 29.97880 PEG 

OflEGfi- 53. 1 5712 DEG LWOOE- 334.0.3604 Dt'G ! M= 32.162'l5 CFG 

>:-■ 3.0169511+03 KM ! V= 4. 956 7;' 35+0 3 KM ! 3- 3. 3327C0 5+0 3 I f' 

:: --6. 6797490+00 KM,'S! VY- 3.8589430+00 KM/S ! VZr 3. 1440.-«- 3-01 Th/S 
PS ELEMFhir'.. 

. 6293S5O4+00..< -. .5890.4845 -npl . 51475927+00? . 5447:..*i59+0C6 

. '..16.64S76+0F5 .. 1 42 + 22 + 0oC 1 0570999 ♦ 00 < . 29V6->- •8 + 0C2 

INPI'T' PhTO IJ-'ING NfiMEI TST 'tINPUl 
■'EOF 

+ + 40PMAI. program TFPMINHtlPP 


OMGINAL vagt . ts 

OF POOR CiUALITY 
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Ia fche 1 iAfe«r*e,6blve bkh^s, the progpem solieits input ttafca with tht ineeeagaej 
a, INPUT MTA USING NAMSL,ISf *$1NP0T' 
b- ENT1H.I C s CQNfINDBi S a STOP 


In bafceh mode, the date input muet be followed by one oard oonfcaliiing the letter* 
as shown. In fche following exantpleei 


Example (1) 
/K^INPUT 


Example it) 

filename. data element 


neoessary data cards 
(see seotion 1 . 1 ) 


ohangee to data element 
(if any) 


.hlEND 

C 

h^LV’UT 


^$END 

C 

gABD filename . data element 


etc. 


eto. 


@E0F gEOF 

SPIN fiFIN 

(h' is a blank that must be included) 


4.P SUBROUTINE PACKAGE 

Along with the ASOP stand-alone program, there is a subroutine pankage that 
may be Included in the u.ser's software. This paokage Is in the form of a 
relocatable element and is located in FM6-N08B69*NUMEG, A30P-5UB. Its opera- 
tion is identical with any user-written subroutine. 

This section will describe the information needed by the user to insure proper 
in.sertlon and operation of the ASOP package within the user's software. 


®The card containing the letter C may bo replaced by just a blank card. 



^•2,1 Suirottfeiiieg 


Tte ASOP aubroiititie paGtege cqnsisbs of seven subf oublnes . These ane a driver 
subroutine (ASQ'P) that aonbrols the basio logie of tbe package and six general 
subroutines that perforin the funations necessary to the analytical theory s 
Th-Sfi© subroutines are 

ASOP Driver subroutine 

CONST Planetary and raatbematical oonsbantS 

P0TJ2 Determines Jg potential 

PSMS PS analytical theory 

PSTOX Transformation subroutinerpS elements 
bo Cartesian coordinates 
TIMEPS Time iteration stopping procedure 

XTOPS Transformation subroutine: Cartesian 

Qoordinates to PS elements 


The subroutines listed above are fully deBoribed in section 3 -Si and a diagram 
of the data flow between these subrautiues can be found in appendix G. 

To help the user add these subroutines to his own software, a relocatable ele- 
ment has been formed that includes all the above subroutines. Therefore, the 
user needs only to include the element 


PM6-N08569»NUMBG . ASOP-SUB 


when forming an executable element, 


4.2.2 Interface Requirement 

To access the ASOP subroutine package , the programer must use the FORTRAN 
iterjent 


CALL ASOP (X,STOP,ISTOP,NEWX) 


A full description of the argument li.'^t variables can be found in sections 
4,2.3 and 3-3 *3- Also, the user must Initialize certain COMMON block variables 
before entering the ASOP subroutine, 

TabLe III gives a list of the variables that must be initialized prior to calling 
the ASOP subroutine and the COMtiON blook in which the variables are located, A 
oompiete description of the var'iables and their allowed values can be found in 
seel Ion 4.1.1. 
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Tms m.~ COMMON BLOCK INITIALIZATION 


Variable 

COMMON 

bloak 

Default 

value(s) 

Page 

lU'NITS 

CPRINT 

1 

61 

CD^ 

DRAG 

0,0 

61 

ABEA® 

DRAG 

0,0 

61 

XMASS^ 

DRAG 

0.0 

61 

ILONO 

PERTRB 

1 

63 

DATE^ 

EPOCH 

5^f0.0 

61 

XJDATE^ 

EPOCH 

0.0 

61 

NMAX^ 

TESS 

2 

611 

MMAX® 

TESS 

0 

6i1 


To be implemented. 

^1,2.3 Input/Outpub Desoripbion 

The argument list bo the ASOP subroutine oonsists of four arguments given 
in the following order. 

CALL ASOP (X,STOP,ISTOP,NEWX) 


On input the variables are 

X An array of eight elements corresponding to the initial state vector 
in the following order: 

X(1): X position 

X(2h I position 

X(3): Z position 

X(^J): X velocity component 

X(5): y velocity component 

X(6): Z velocity component 

X(7): Physical time (set to aero) 

X(8): Total energy (set to zero) 
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STOP stop value desired] it may be given in days or ravolufeions. 

ISTOP Flag determ i n i ng wtvetber the value given to STOP is in days CISTOP 
= I) ©r revolutions (ISTOP ^ , 

NSWX Flag determining if the A50P suhroufeine is to be initialized 
NO ^ 0, ISS B 1 


The initialization process ®asb be done whenever new initial 
Qonditions are entered. 


WARNINCi! A npraber must not he used in this position of the argument list, 

Assign the desired value to a variable, and use the variable in the 
argument list, if a number is used in this position instead of a 
variable, unpredictable results may occur. 


Input to the ASOP program is also done by means of COMMON block's, These COMMOW 
block input variables control the internal operation of the ASOP subroutine 
paokage and should not be changed once the subroutine has been initialised. 

Table III (seation ^1.2.2) gives a oomplete list of the variables that must be 
initialized and gives their default values. 

On output the variables are 

X An array of eight elements oorresponding to the final state vector 
at the given value of STOP. The order is the same as for input. If 
the value of STOP was given in days (ISTOP = 1), then the value of 
X(7) will be an approximate value of STOP given in the time units 
specified by lUNITS. 

STOP ) 

> Onohanged from the input values, 

ISTOP ) 

NEWX Set to zoi-o (0) upon return from the ASOP subroutine; therefore, 
a variable name should always occupy this position in the argu- 
ment list. Otherwise, unpredictable results may ooour. 
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APFBH'DJX. A 


A7AILABLB ONiyS ANB PHYSICAl, COHSTAMTS 

WifchiD tibe A$QP progpunv, tbere are a numtoep of* opbiona for* bhe itipufe upibs, 
The ooiapabibility of fcbo Input values with the selected physical constants 
is the respoaeiblllty of the user and is controlled by bhe Input flag IDHITS. 
The followingvQOnstaabs are listed in order according to the value assigned 
bo IttNlTS. If no value for lUNITS is explicitly given, then 1 is assumed. 


= Earth radius (equatorial) 


= 6378.140 

km 

(1 NITS = 1,5) 

= 3443.920 

ntn 

(lUHITS = 2,6) 

= 2,092566 

X 10^ ft 

(lONlTS = 3) 

= 6.378140 

X 10^ m 

(IDNITS = 4) 

= 1.0 E.r. 


(lONITS = 7) 


U = Gravitational constant of the Earth 


= 3*986013 X 
= 6.ET5029 X 
= 1.4076*17 X 
= 3-986013 X 
= 5.165873 X 
= 8.132438 X 
a 5.530432 X 


10^ kra^/seo^ 

4 3 2 

10 run /sec 

lo’'® fb^/seo^ 

10^^* in^/seo^ 

10^^ kra^/hr^ 

10 ” nra^/hr^ 

10"3 E.r.3/tnin^ 


CIUNITS - I) 
(IWJTS = 2) 
(la’Hifii = 3) 
(lutm's = 4) 
(TUNITS - 5) 
(lUNITS = 6) 
(lUNITS = 7) 


DA7 = Time conversion 


10^ 

sec/ day 

(IONIT3 = 1,2, 3, 4) 

10^ 

rain/ day 

(lUNITS = 7) 

10^ 

hr/ day 

(lUNITS =5,6) 
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J 


s 


e 


Jg ooeffioleofc or the Earth's gtopobenbial 
1.082637 X 10“^ 


E EPS a 3/2 (y Jg R/) 

Tf = PI = 3.14159 26535 898 


If Keplerian e3,en:enfcs are Ihput (lEL = 1), then the semiraajor axis must be 
iapufc in a distance compatible with the selected value of lUNITS, i.e., if 
ItJHITS s 2, then EL(1) must be given in nautioal miles; if lOtJITS = 4, then 
EL(1) must be In meters, etc. All angles must be given in degrees, 

For Cartesian coordinates, the input unit must be 


IDNITS EL(1) ^ EL(3) EL.(4) ^ EL(6) 


1 

kra 

kra/ sec 

2 

nm 

nm/seo 

3 

ft 

ft/seo 

4 

m 

m/seo 

5 

km 

km/ hr 

6 

nra 

nm/hr 

7 

E.r 

E.r/min 


All computations within ASOP are done using the input units. 
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Appmeg B 


BiQaSEBD CARBS 

1 . fM6,«H09S69 

a, >|ftS&jA *NtlME)G. 


3 - 

4 . 

5. 


6 . 


>iXQT »WMEG.ftSOP^PRGG 
INiPGT M'SA USING NAMELIST '$INPUT» 
»^$INPUT . 

Keoessary input paramebens; [ 


see seetion 4.1.1. All ) on { 


parauiebers must be pneoeded 
by at least one spade* 

/ 


>gADD filename.elenient name 

> Input the neaessany oonreotion 
if needed. 


Initial output 

7, ENTER: C = CONTINUE j. S = ST0P 


If an S is entened in response bo the prompt then the message 

»»NORMAL PROGRAM TERMINATION*# should appear. 

If a C is entered, the program will printout the desired information, 
and the sequenoe will begin again at line 

8. >@E0F 

**N0RMAL PROGRAM TERMINATION** 

> = system prompt; 
if = neoessary space 


ORIGINAL J’AGB B 
OF FOQR (^UALlTy 
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Agpp mnuLT uiOEs 

lEL = 1 (Keplerian element 3 ) 

EL = tJone j must be input by user 

STOP = 100.0D-0 

ISTOP = 2 (STOP value is in revolutions) 

PBINT = O.ODO 

IPRIHT = 0 (No printout; PRINT value ignored) 
IDEAG^ - 0 (No drag terras desired) 

AREA^ = O.ODO 

CD® = O.ODO 

XMASS^ = O.ODO 

ILONG =1 and sho period secular terms) 

NMAX® = 2 (Include J sonal terra) 

|| u 

MMAX^ =0 (No tesseral terms) 

lELPRT = 0 (Do not print PS elements) 

lUNITS = 1 (Input and output values are given as 

distance = kra 
velocity = km/seo 
time = day 
angles = deg) 

®To be implemented. 
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AFPEWDIX D 


aoiEOQTItjTB gYOLB TXMBS Md STORAGE RjiaBiaMBKI'S 

AU ttijieB refer te fche eaceeufciofl of a fORfMN ? profraw oo a MIllAfi 1 1 ia*^iX.BS 
eystem:, 

fOBTRAW oompileri S0E3 
MAP prooetaor t gy.! RL71-3 

gborafie reqT 4 Apen)eo.bs , words 


Subroutine 

name 

Cycle 
timear ms 

Octal 

Deoimal 

MAIN 

®34 to 42 

145 

101 

AEIXY 2 

<1 

356 

238 

ASOP*^ 

2 to 11 

134 

92 

CONST^ 

<1 

207 

135 

INPUT 

— 

272 

186 

OUTPUT 

aS to 33 

50 ? 

326 : 

P 0 TJ 2 *^ 

<1 

4 l 

33 

PSANS*^ 

1 

1773 

1019 

PSTOX^ 

<1 

360 

240 

TIMEPS^ 

4 to 9 

202 

130 

XTOPS*^ 

<1 

540 

352 

XYZAEI 

<1 

406 

262 

Totals 

A 30 P program 

ASOP subroutine package 

6047 

3721 

3111 

20 f '1 


S 

Does not inolude the time required for i \ Input. 

“ASOP subroutins package programs. 

OF POOR 



The oyole tiroes showe are the total times needed >y the module to return an 
answer after it is oalled. Therefore, if a module roust oall a suhroutiDe, 
then the time shown for that module la the exeoution time needed by the module 
and the subroutine nailed- 


Bxaroples The time of 2 to 11 railliseoonds given for the ASOP suhroutine 
is the time required by the module to 

a. Aooept the input data 

b. Transform the input data into the ps elements (oall XTOPS) 

0 , Initialize the analytical theory and set the necessary constants 
(call PSANS and G0N9T) 

d. Compute the new PS elements (oall PSANS) 

e. Transform the new PS elements baok into the Cartesian coordinates 
(oall PSTOX) 

f. Return to the calling routine 


Values given for the storage requirements of the individual subroutines are 
values returned by the FORTRA.N compiler when forming a relocatable element. 
The final, executable program will require more space because of the system 
library modules that must also be inoluded . 
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APPBMPIX S 


gEMBPAL VARlAaLB ABBRSVIATIOaS ASP DBPIMITIQWg 
eletneQt>a i 

Coord ingfc$f cr^» cr^i 0j> Oj| 

(KobSJ 0g ~ p^, C Pg, CtJ - P^j 0g =* P|j) 

Momenba P-ji pgj Pg* Pi} 

IndepeMepb variable T (true anomaly) 


PS Hatnilfconiam 


P 



r2 


V 




Where q = G - V2 $ + — — ■, and V is the perturbing potential. 

2 2L 


PS elements : 

(j) = true anomaly 
g = argument of perigee 
h = argument of the asoending mods 
2, = time element 

$ a conjugate to «(i, related to the two*-body energy 

0 a total angular momentum 

H a Z component of the angular momentum 

I, a total energy (two body plus perturbing potential) 

(Note; For a oomplete deacri ptioo of the relationship between the DS and 
the PS elements, see reference 2.) 


OmGIN AL P 
Of POOE QUAEl-'J^' 
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't “ Xy, Xg) 3 pggitiQii veolor 
a (Vj^j Yy, Vg) = veXooifey veotop 
r s magnifeude of fehe pasibiOB veofcop 
b a physioal time 

KdPlerlaH elaaeptg j 

a a semimajop axis 

e = eedentPioity 

i a inolioabian to the equator 

y = argunieab of perioenbep 

j2 = argument of the ascending node 

M a mean ano®aly 

Planetary variables : {see appendix A for the numerical values used) 

Hg = equatorial radius 
y = gravitational oonstant 

General : 

kra s kilometers 
nm = nautical miles 
ft = feet 
m = meters 
E.r. = Earth radius 
sec a seconds 
hr a hours 


rain = minutes 
rad = radians 
deg = degrees 
t = time 

^ = denotes a vector as X 


90 



APPENDIX P 


BQDA'flOltS or THE ANAqTIGAL TRECRX 

A oomplefce first order solution for the motion of a setalllta perturbed by o''- 
lateness has been developed (ref. 2). A brief outline waa given in refaranoe 
6 end is reproduoed in this appendix. 

The Hamiltonian for the J* perturbed oasa oan be written as 


U 

F e p. - -jss + e F 

where 



s = 3/2 P r/ 


R is the mean equatorial radius of the oentral body; y Is the gravitational 

constant of the central body, and J_ ia the oblatenesa ooeffioient. 

d d 

The differential equations are solved by a method of Von-Zeipnl. The elements 
undergo a oanonioal transformation through a determining function so that 

the short periodic terms are eliminated from the Hamiltonian. The equations 
of motion in the transformed system may then be solved with an accuracy 
of order c. 

The solution algorithm can be divided into three steps: 
a. Canonical transformation to eliminate the short periodic terms: 


33 1 


'’k,0 °k,0 ^ 


9Pl<,0 




3Si 

^’>^.0 ' '" x.o “ '" o ’ 




ORIGINAL PAGE IS 
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b, The analytical Intsgmblen of the fenanifopwed squatien* of BOtlatj; 

Cff c cri ^P + A^t 

Q oos CAgt) ein (AgT) 

"^3 " tAjT) 

■ '^4,0 + 

"Pl,Q 

n = 4,0 •*’ ^2,0 

p^ = ooe (AgT) + aj^Q sin (AgT) 

n "4,0 

o. The baok transf ormafcloti j 

9S^ 

k * 1,2, 3,1 

If one defines 

3Si 
3crj^ 

than 

1 

®ik = 5 “k j 

whi-.'re 

Q 

w = — - 
2pq 


2 w y Gi^ 

0 ■■' 
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1 






3 






(fj/ * 



= 0 far If = 1,3, 4,7*3 

Qg =s -0g 



1 


h = ~^,6 

Pr Pkv <1> qjt* n^’ %’ displayed 


1 


r 1 1 

2 h\ 

4 JL. ' 

XI 

II 

TT \ 


i 

■^^ 6 ] 






2-M^ 

C2ag)3^2 


p^ = 0 foP k = 1,3,43,7 


1 2 2 
q a - - (0g + rig 



am^^-^fXB 

of POOB QpAXdl^ 
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<?« “ 




5 " i 


^8 " 


2 (affg)3/2 

q,{^ ™ 0 tor k= 1,3, 4, 6, 7 


Q = 


71 


i - - (sf + o^h 


1/2 


*^2 ^ 

2QM2 


^8^6 


Qo = 


VW 


= 0 for k ^ 1,3, 4, 3, 7 


S, r 


5 Og - 1 (Oj,! - U^g) 


B “6 , 

*12 = i * 5- - i t'’6°2 - Va^ 


D c 


^6 


1 


= “+ Bfi“ - (CTcOr - SJoHr) 

Id 3 2 3 “ 2 0*5 2f> 


1 

'^ik = T~ \ ~ o ^‘^6°k - ^ = '>,3,4,5, 7, 8 

J w 


’'1 ■ 5 “a * 5 <V - V' 
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/ 1 e\ 




i/ 3 


1 


■* ? + J?g ^? g 

IS 2 I ® 2 


J ^5 ^ 

■ Bg + a 


■r.,. = 2 s t 1 (sjs^ + Bgo^ 


Ik - 3 "k " a 


a) for k e 


s. = - 


1 /o 


^Sk " 2Q t| % 


for k “ 1 )2,3 1 . ‘ • i' 


Y2 ^ m 


Y - - 1,^(- Qr. - Oi \ ^or k ^ l,2*3j* ' *!® 

%k “ aqU ^ 


'J 

§2 = - g (ffgS 0$G) 


«sa = - J 'Vs ♦ V2 * 

636 ^ “ I ^"2®6 + 

«ai< = - i '“J'k * “6»k' ‘■°’' 

7j = I - "j") 

■^32 ' I ''’6“2 - '’a“a ■ “’ 


k s i,3.i».5i7.8 


OETGmL PAGE IS 
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’^36 " i '“6“6 - 

■^3k = I <"8^ - Vk> ‘‘•=^ ' 

Tij^ ^ sin ^,0 ^ 

:; cas Aa^ for JL - 1>S»3 
Hare c, 6, Oj^, B, H and H, 


0 = (G + H) 


0^2 - 02^ 


H30 


3= t^- “=*» "3 


HyC 

c = -■ - ■ + CG + H) cr 

7 (Q + H) " 


+ \ 

C, = -7~~ c for k = 1,2, 

k (G + H) 


s = -(G + H) CJ3 dy 
H3S 


3„ ~-r. ^ - (G + H) 0n 

3 G + H ' 


HyS 

3 = CG + H) 0 

7 G + H 3 


(G{^ + 

a 3 ^ 3 for k s 1 j 

k (G + H> 

B = G^ - 3H^' 


Bj^ = 2CGG|^ - 3 HHj^) for k 


= 1,3, '1, 5, 7, 8 


are di»iDlayed 


4, 5, 6, 8 


2,*1,3»6,8 


1,2,3 ,. ..,8 
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W s § - I (03^ + 


H3 S -^03 


s -I 


7 = -“7 


®k ■ °l£ '' ° Ii2i'ti5j6,8 


AtsbreTiafeioBS used in the Integrabion th« prijnsd ayntein 


- 3’ * 


‘3 = P3 


= I 


I b - ^ + fbg 


+ 


= 1 + 


(t- |)+A. 


r a - — 

pq 


i’l = QuP + 2q 


r s (up + 

^ u 


' { 2 P„) 3 / 2 V 2 


I'l, * 


Up <■ ^vw 



.2 


b = 1 
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MAIN 



OP POOB QUALITY 


hasa-jsc 













